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This  is  the  final  report  on  the  work  done  in  the  area  of  high  resolution 
photoacoustic  and  photothermal  imaging.  It  contains  recent  advances  in  photoacous¬ 
tic  and  photothermal  theory  and  the  experimental  demonstration  of  new  techniques. 
Photoacoustic  and  photothermal  theories  have  been  extended  to  include  the  effects 
of  the  highly  focused  optical  power  sources  necessary  for  high  resolution  imaging 
and  three  high  frequency  techniques  (1  GHz)  have  been  demonstrated  and  used  to 
characterize  the  material  properties  of  solids. 

The  formalism  behind  photothermal  characterization  of  solids  is  established 
under  general  focusing  conditions.  >An  analysis  of  the  frequency  dependence  of  the 
photothermal  response  demonstrates  that  very  high  frequencies  may  be  employed 
without  loss  of  sensitivity  if  highly  focused  optical  beams  are  used.  It  is  shown  that 
the  photothermal  signal  can  be  used  to  directly  measure  the  thermal  conductivity 
of  solids  on  a  microscopic  scale. 

The  three  dimensional  thermoacoustic  field  equation  is  derived  and  discussed.  > 

/ 

This  equation  is  used  to  obtain  the  one  dimensional  thermoacoustic  wave  equation. 
It  is  shown  that  an  error  exists  in  the  traditional  one  dimensional  equation  used 
in  photoacoustic  theory.  A  new  approach  to  the  solution  of  the  one  dimensional 
problem  is  presented.  It  is  shown  that  the  far  field  stress  can  be  related  to  the 
thermoacoustic  source  distribution  by  a  Fourier  transform  relation  under  arbitrary 
acoustic  boundary  conditions. 

Two  high  resolution  photoacoustic  techniques  have  been  demonstrated.  These 
techniques  involve  the  use  of  an  acoustic  lens  for  collection  of  the  acoustic  power 
generated  by  a  highly  focused  optical  beam  modulated  at  1  GHz  frequency.  Images 
of  gold  and  laser  recrystallized  silicon  films  are  presented  demonstrating  the  high 
resolution  and  sensitivity  of  these  techniques. 

A  high  resolution  photothermal  probe  has  been  demonstrated.  This  optical 
technique  provides  a  means  of  photothermal  characterization  with  sub-micrometer 
resolution  and  high  sensitivity.  The  theory  behind  the  probe  is  presented  along 
with  the  experimental  verification.  Images  of  boron  implanted  silicon  are  included 
which  demonstrate  the  resolution  and  sensitivity  of  this  technique. 
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Chapter  1. 


INTRODUCTION 

Over  the  past  three  decades,  a  substantial  interest  in  the  microscopic  properties 
of  materials  and  material  systems  has  been  demonstrated.  Many  new  techniques 
and  tools  have  been  developed  for  microscopic  examination.  The  use  of  the  scanning 
electron  microscope  for  elemental  analysis  on  a  microscopic  scale  (microprobe)  has 
provided  useful  information  about  the  composition  of  very  small  systems  [1.1).  Work 
done  in  the  area  of  transmission  electron  microscopy  has  brought  the  state  of  the 
art  to  the  point  where  individual  atoms  can  be  resolved  [1.2].  Research  in  the  area 
of  acoustic  microscopy  has  extended  the  resolution  of  the  water  coupled  instrument 
beyond  that  of  optical  microscopy  [1.3].  In  cryogenic  liquids,  the  highest  resolution 
demonstrated  to  date  is  approximately  200  angstroms  [1.4].  There  is  hope  that 
further  improvement  in  resolution  may  be  possible.  Scanning  optical  microscopes 
have  been  pushed  toward  their  limits  using  short  wavelength  laser  sources  to  obtain 
the  highest  resolution  [1.5].  Recently,  a  near  field  approach  has  extended  the  the 
resolution  of  optical  microscopy  far  beyond  the  diffraction  limit  [1.6].  Tunneling 
probes  have  been  successful  in  imaging  the  electronic  properties  of  materials  with 
resolution  on  the  atomic  scale  [1.7].  All  of  these  efforts  are  manifestations  of  the 
interest  in  the  microscopic  properties  of  materials. 

Two  areas  of  microscopic  imaging  which  have  received  attention  in  recent  years 
are  photoacoustic  and  photothernial  imaging  [1.8-1.14].  Photothermal  imaging  is 
defined  to  include  those  techniques  which  are  capable  of  independently  measuring 
the  thermal  response  of  a  solid  generated  by  absorption  of  optical  power.  The 
information  obtained  by  these  methods  does  not  contain  acoustic  artefacts  which 
may  complicate  the  interpretation  of  the  data.  Photoacoustic  techniques,  on  the 
other  hand,  are  those  which  are  inherently  tied  to  the  coupling  of  thermal  and 
acoustic  waves  in  the  material  to  be  imaged.  These  techniques  can  be  used  to 


measure  both  the  thermal  and  acoustic  properties  of  materials.  Each  of  these  areas 
has  been  usefully  applied  to  microscopic  examination  of  materials  and  systems. 

Photoacoustic  and  photothermal  imaging  are  based  upon  the  generation  of  a 
periodic  thermal  distribution  in  solids  by  absorption  of  modulated  optical  power. 
See  Figure  1-1.  The  thermal  response  can  be  treated  as  a  sum  of  thermal  waves 
generated  by  a  distribution  of  point  sources  in  the  region  of  optical  absorption. 
These  thermal  waves  are  temperature  waves  whose  characteristics  depend  upon  the 
thermal  properties  of  the  medium  in  which  they  are  generated.  Using  highly  focused 
optical  beams,  this  dependence  can  be  exploited  to  obtain  information  about  the 
thermal  properties  of  a  solid  with  microscopic  resolution. 

There  are  several  reasons  for  the  interest  in  the  microscopic  thermal  properties 
of  materials.  In  crystalline  semiconductors,  the  thermal  conductivity  is  strongly 
dependent  upon  the  quality  of  the  crystalline  lattice.  Perturbations  to  the  crystalline 
structure  caused  by  defects  or  dopants  can  greatly  modify  the  thermal  conductivity 
of  the  material.  The  measurement  of  this  thermal  parameter  is  therefore  a 
means  of  characterizing  the  crystalline  structure  as  well.  Another  source  of 
interest  is  in  the  area  of  subsurface  characterization.  The  thermal  waves  can 
provide  useful  information  about  subsurface  structure  such  as  cracks,  delaminations 
and  interfaces  that  cannot  be  evaluated  by  other  non-destructive  techniques. 
The  ability  to  microscopically  measure  the  thermal  properties  of  materials  and 
non-destructively  evaluate  subsurface  properties  provides  useful  information  for 
materials  characterization. 

Because  thermal  waves  are  highly  damped,  they  are  difficult  to  detect  directly. 
Therefore,  various  indirect  ways  for  probing  the  thermal  response  have  been 
developed.  A  diagram  of  the  generic  classes  of  these  techniques  is  shown  in  Figure  1- 
2.  As  can  be  seen  from  the  figure,  the  techniques  can  be  classified  into  three  general 
classes,  thermo-optic,  thermoacoustic,  and  pyroelectric.  Each  of  these  techniques 
has  its  advantages  and  disadvantages,  depending  upon  the  sensitivity,  simplicity, 
and  spacial  resolution  required.  Most  of  the  techniques  are  generally  useful  at 
modulation  frequencies  below  one  hundred  kilohertz,  where  the  thermal  wavelength 
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.  Many  of  the  available  thermal  wave  imaging  techniques  divided 
into  generic  classes. 


in  a  semiconductor  is  greater  than  10  microns.  Some  experiments  have  been  done 
at  10  MIIz  frequencies  where  the  thermal  wavelength  in  a  semiconductor  is  greater 
than  1  micron.  The  work  which  will  be  presented  here  has  been  done  at  frequencies 
near  1  GHz.  The  thermal  wavelength  in  these  experiments  is  somew-hat  greater 
than  1000  angstroms.  Under  these  conditions,  the  resolution  of  the  thermal  imaging 
system  is  essentially  limited  by  the  optical  spot  size. 

In  order  to  understand  and  evaluate  the  photoacoustic  and  photothermai 
experiments  that  will  be  presented  in  the  last  two  chapters  of  this  work,  Chapter  2 
has  been  devoted  to  understanding  the  photothermai  response  of  materials.  This 
requires  an  introduction  to  thermal  waves,  and  an  analysis  of  how  the  thermal, 
optical  and  geometrical  properties  of  materials  are  related  to  the  photothermai 
response.  A  discussion  of  the  applications  where  photothermai  characterization 
may  be  useful  will  then  be  presented. 

Chapter  3  contains  an  analysis  of  one  and  three  dimensional  thermoacoustic 
generation.  This  theory  is  useful  for  understanding  and  comparing  the  various 
thermoacoustic  techniques,  two  of  which  will  be  presented  in  Chapter  4.  The  theory 
will  show  that  the  thermoacoustic  conversion  efficiency  increases  with  modulation 
frequency  for  thermal  sources  which  are  small  compared  to  an  acoustic  wavelength. 
Also,  it  will  be  shown  that  the  existing  theory  of  one  dimensional  thermoacoustic 
generation  predicts  a  generation  efficiency  which  is  nine  times  smaller  than  that 
predicted  by  the  correct  theory. 

In  Chapter  4,  two  experiments  will  be  described  which  involve  an  acoustic  lens 
as  a  collector  of  photoacoustically  generated  sound.  The  merits  of  using  a  focused 
acoustic  collector  will  be  discussed  and  the  sensitivity  of  the  acoustic  lens  as  a 
collector  of  photoacoustically  generated  sound  will  be  calculated.  Images  will  then 
be  presented  demonstrating  the  sensitivity  and  resolution  of  these  techniques. 

Chapter  5  contains  the  description  of  a  new  photothermai  laser  probe. 
The  theory  behind  the  probe  will  be  presented  and  a  theoretical  sensitivity 
calculated  under  reasonable  experimental  conditions.  A  figure  of  merit  will 
be  established  for  evaluation  of  various  acousto-optic  coupling  media.  The 


probe  will  then  be  compared  with  two  other  high  resolution  techniques.  The 
experimental  implementation  will  be  described  and  verification  of  the  theory  will  be 
reported.  Finally,  photothermal  images  of  boron  implanted  silicon  will  be  presented 
demonstrating  the  sensitivity  and  resolution  of  this  technique. 
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Chapter  2. 


PHOTOTHERMAL  RESPONSE  OF  MATERIALS 


2.1  Introduction 

Both  photoacoustic  and  photothermai  imaging  have  their  basis  in  the  conver¬ 
sion  of  modulated  optical  power  into  a  periodic  temperature  distribution  in  a  solid. 
The  thermal  response  depends  upon  several  parameters  which  include  the  charac¬ 
teristics  of  the  heating  beam  i.e.  modulation  frequency  and  spot  diameter,  and 
the  material  properties  and  geometry  of  the  absorbing  system  i.e.  optical  absorp¬ 
tion  constant,  thermal  conductivity,  density,  specific  heat  and  film  thickness.  An 
understanding  of  these  factors  makes  possible  the  characterization  of  solids  by  ex¬ 
perimental  measurement  of  the  photothermai  response. 

In  this  chapter,  the  photothermai  response  of  solids  will  be  examined.  Since 
this  response  can  be  treated  most  clearly  using  a  thermal  wave  (or  Green’s  function) 
approach,  a  discussion  of  thermal  waves  will  be  presented.  The  dependence  of 
the  thermal  response  upon  the  modulation  frequency,  the  optical  and  thermal 
properties  of  the  solid,  and  geometry  will  then  be  examined.  In  the  final  section, 
the  possibilities  for  useful  application  of  these  principles  will  be  discussed. 

2.2  Thermal  Waves 


The  thermal  response  of  a  system  is  governed  by  the  inhomogeneous  thermal 
diffusion  equation  shown  below  (2.1). 


V20- 


IdO 
k  dt 


(2.1) 


The  variable  0  is  the  temperature,  k  is  the  thermal  diffusivity,  K  is  the  thermal 
conductivity,  and  A{?,t)  is  the  power  per  volume  released  in  the  medium.  If  the 


source  function  A(r,  f)  is  harmonic, 


A(?,t)  =  A(r)e}ujl  (2.2) 

the  diffusion  equation  becomes 

V20  +  k26=^A(r),  (2.3) 

where 

*•  =  (£)**' +J>-  <2-4» 

This  equation  looks  like  the  Helmholtz  equation,  encountered  in  electromagnetics. 
However,  the  term  kf  is  complex  in  the  thermal  case,  whereas  in  the  electromagnetic 
case  it  is  real. 

The  Green’s  function  for  this  equation  can  be  obtained  if  we  substitute  a  delta 
function  for  the  spacial  dependence  of  the  source  function. 

A(r,t)  =  6(r-?0)e>’“t  (2.5) 

With  this  as  the  source  term,  the  solution  is  given  by  [2.2] 

=  >2-6' 

This  is  the  Green’s  function  solution  to  the  thermal  diffusion  equation  with  a 
harmonic  source.  This  solution  has  been  called  a  thermal  wave.  One  of  the 
advantages  of  a  harmonically  varying  thermal  source  when  compared  to  a  pulsed 
source  is  that  thermal  waves  have  a  fixed  velocity  and  wavelength.  Under  these 
circumstances,  a  phase  delay  can  be  related  to  a  distance  or  a  variation  in  the 
thermal  properties  of  the  medium.  In  contrast,  when  a  pulse  of  heat  is  generated, 
all  the  Fourier  components  of  the  heat  pulse  propagate  at  different  velocities,  and 
the  pulse  shape  changes  rapidly  with  the  diffusion  of  heat  to  the  surrounding  areas. 
Therefore,  no  unique  velocity  can  be  attributed  to  the  transfer  of  heat. 

The  periodic  thermal  response  at  any  point  due  to  aspacially  distributed  source 
can  be  obtained  by  integrating  the  product  of  the  source  function  and  the  Green’s 


function  over  the  volume  where  the  source  is  finite. 


Ge(r\?0-u)A{r0)dV0 


(2.7) 


There  are  several  interesting  properties  of  thermal  waves.  First,  these  waves 
have  k-vectors  which  have  equal  real  and  imaginary  parts.  They  are  therefore  very 
highly  damped.  Second,  the  thermal  k-vector  is  proportional  to  the  square  root  of 
the  frequency.  By  contrast,  in  the  frequency  range  of  interest,  the  acoustic  k-vector 
is  proportional  to  the  frequency.  This  results  in  a  change  in  the  ratio  of  the  thermal 
to  acoustic  wavelength  as  the  frequency  is  changed.  At  low  frequencies,  the  thermal 
wavelength  is  very  small  compared  to  the  acoustic  wavelength.  At  high  frequencies, 
the  thermal  wavelength  begins  to  approach  the  wavelength  of  sound.  However,  at 
frequencies  where  the  thermal  wavelength  approaches  the  acoustic  wavelength,  the 
thermal  waves  begin  to  lose  their  diffusive  character.  This  is  true  because  equal 
wavelength  implies  equal  velocity.  If  the  velocity  of  the  thermal  wave  approaches 
that  of  sound,  it  cannot  be  diffusive.  Fortunately,  for  most  materials,  this  occurs 
at  frequencies  near  100  GHz,  Far  From  the  region  of  interest  (1  GHz). 


Third,  materials  can  be  characterized  by  a  thermal  impedance.  The  thermal 
impedance  Z,-  is  given  by  the  ratio  of  the  temperature  to  the  heat  flux  for  a  traveling 
thermal  plane  wave  in  medium  ».  This  thermal  impedance  is  related  to  the  material 
properties  by  the  following  relation. 


(2.10) 


See  Figure  2-1.  Using  this  impedance,  a  thermal  reflection  and  transmission 
coefficient  can  be  obtained  using  the  following  relation. 


(Z2-Z,) 
(Z2  +  Z\  ) 


(2.11) 
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Figure  2-1.  Illustration  of  the  reflection  and  transmission  of  a  thermal  plane 
wave  that  takes  place  at  an  interface  between  two  media  of  Afferent 
thermal  impedance.  The  thermal  impedance  of  medium  i  depends 
upon  the  radian  frequency  *>,  the  thermal  conductivity  h;,  the 
density  and  the  specific  heat  C,\ 


Because  the  generation  and  propagation  of  thermal  waves  depends  upon  the 
geometrical,  thermal  and  optical  properties  of  the  medium  in  which  they  exist, 
a  measurement  of  the  thermal  response  makes  possible  characterization  of  the 
medium.  In  the  next  section,  the  dependence  of  the  thermal  response  upon  several 
parameters  will  be  examined.  These  results  demonstrate  that  information  about 
material  properties  can  be  extracted  from  the  photothermal  response. 

2.3  Characterization  of  Solids  by  Photothermal  Response 

The  photothermal  response  of  a  solid  depends  upon  many  different  geometric 
and  material  parameters.  The  important  geometric  parameters  are  the  diameter 
of  the  optical  spot,  the  thickness  and  geometry  of  the  materials,  and  the  ratio 
of  these  parameters  to  the  thermal  wavelength.  The  material  parameters  which 
effect  the  photothermal  response  are  the  thermal  conductivity,  specific  heat,  optical 
absorption  coefficient  and  the  thermal  impedance.  In  order  to  extract  useful 
information  about  a  solid  from  a  photothermal  measurement,  the  dependence  of  the 
thermal  response  upon  these  parameters  must  be  well  understood.  In  this  section, 
these  parameters  and  their  effect  on  the  thermal  response  will  be  discussed. 

When  the  optical  spot  size  is  very  large  compared  with  a  thermal  wavelength, 
the  photothermal  response  can  be  treated  one  dimensionally.  The  results  from 
the  one  dimensional  theory  will  be  summarized  in  the  next  section.  When  the 
optical  spot  size  is  smaller  than  or  on  the  order  of  the  thermal  wavelength,  a 
three  dimensional  theory  must  be  used.  The  three  dimensional  analysis  has  been 
done  using  the  numerical  routine  found  in  the  Appendix.  The  results  of  these 
calculations  are  contained  in  the  sections  which  follow.  Included  in  this  analysis  is 
the  dependence  of  the  thermal  response  upon  frequency,  optical  absorption  depth, 
optical  spot  size,  and  thermal  conductivity.  Finally,  a  discussion  of  the  application 
of  this  characterization  capability  will  be  presented. 


2.3.1  One  Dimensional  Photothermal  Response 


When  the  optical  spot  size  is  many  thermal  wavelengths  in  diameter,  the  ac 
temperature  at  the  surface  of  the  solid  due  to  the  absorption  of  the  modulated 
light  can  be  treated  with  a  one  dimensional  theory.  Most  of  this  theory  has  been 
extensively  described  by  Rosensweig  and  Gersho  [2.3]  and  therefore  will  only  be 
summarized  here.  There  are  five  cases  of  interest.  The  model  used  here  consists  of 
three  regions:  a  vacuum,  an  absorbing  layer,  and  an  infinite  backing.  Modulated 
light,  incident  through  the  vacuum,  periodically  heats  the  absorbing  layer.  This 
region  is  characterized  by  a  thermal  conductivity  K ,  density  p,  specific  heat  C, 
optical  absorption  coefficient  ,3,  and  thickness  L.  The  infinite  backing  is  charcterized 
by  a  thermal  conductivity,  density  and  specific  heat.  The  resulting  harmonic 
temperature  response  at  the  surface  of  the  absorbing  layer  is  calculated  for  an 
optical  intensity  modulation  /<>  at  a  radian  frequency  u>.  The  thermal  k-vector 
is  given  by  kt .  When  the  thickness  of  the  absorbing  layer  is  much  less  than  the 
thermal  wavelength,  it  is  called  thermally  thin.  If  its  thickness  is  much  greater 
than  the  thermal  wavelength,  it  is  called  thermally  thick.  Likewise  with  the  optical 
absorption  depth. 

Case l  -  Thermally  thick,  optically  thin 

In  this  case,  the  magnitude  of  the  sinusoidal  temperature  rise  at  the  surface  of  the 
absorbing  layer  is  proportional  to  the  optical  attenuation  constant,  and  inversely 
proportional  to  the  modulation  frequency,  density  and  the  heat  capacity. 


e=  ~^r° 
2u>pC 


(2-12) 


Case‘2  -  Thermally  thick,  optically  thick,  for  0  much  greater  than  kt 
The  temperture  in  this  case  is  inversely  proportional  to  the  square  root  of  the 
thermal  conductivity,  density,  heat  capacity,  and  frequency.  Note  that  this  is  the 
case  treated  by  White  [2.4]. 
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CaseS  -  Thermally  thick,  optically  thick,  for  kr  much  greater  than  /? 

The  temperature  in  this  case  is  proportial  to  the  optical  absorption  constant,  and 
inversely  proportional  to  the  frequency,  density  and  heat  capactity. 


2uipC 


(2.14) 


Case-i  -  Thermally  thin,  optically  thick,  for  /?  much  greater  than  kT 
In  this  case,  the  temperature  is  inversely  proportional  to  the  square  root  of  the 
thermal  conductivity,  density,  heat  capacity  and  frequency.  There  is  a  complicated 
dependence  on  the  thermal  wave  reflections  at  the  interface  betw-een  the  absorbing 
layer  and  the  substrate. 


— /p(  1  —  j)0  fe‘’L  +  re-oL\ 
2s/2^/K~pCu  +  re~°L  ) 


(2.15) 


<7  =  (S)2(1  +  ;)  (216) 

The  thermal  wave  reflection  coefficient  at  the  interface  is  given  by  r. 

C'gseS  -  Optically  transparent  film  on  an  absorbing  substrate 
If  the  harmonic  temperature  rise  at  the  interface  between  the  transparent  film  and 
the  absorbing  substrate  is  given  by  0,-,  the  harmonic  temperature  at  the  surface  of 
the  film  is  proportional  to  an  exponential  with  an  exponent  which  is  the  product  of 
the  thermal  k-veetor  and  the  thickness  of  the  film.  See  Figure  2-2. 


0  =  9ie~k,L  (2.17) 

There  are  several  interesting  points  that  can  be  made  with  the  results  presented 
above.  First,  it  is  apparent  that  if  one  desires  a  sensitivity  to  the  thermal 
conductivity  of  a  material,  the  optical  absorption  depth  should  be  short  compared 
to  the  thermal  wavelength  (/?  »  kt).  This  seems  intuitive  since  no  thermal  diffusion 
takes  place  in  a  region  which  is  uniformly  heated.  Second,  the  product  of  density 
and  specific  heat  seen  in  many  of  the  expressions  above  does  not  contain  very 
much  useful  information,  because  at  room  temperature  this  product  varies  little 
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Figure  2-2.  Illustration  of  how  the  surface  temperature  of  a  film  depends  upon 
its  thickness  /,  and  the  thermal  k-vector  kt  when  the  substrate  upon 
which  is  sits  is  harmonically  heated. 
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from  material  to  material.  Third,  Case  4  indicates  that  the  thermal  properties 
of  a  subsurface  boundary  between  two  materials  may  be  evaluated  by  measuring 
the  reflected  thermal  waves  from  the  interface.  High  frequency  thermal  waves 
are  uniquely  suited  to  such  subsurface  characterization.  Fourth,  Case  5  indicates 
that  the  measurement  of  film  thickness  may  be  possible  when  the  film  is  optically 
transparent.  Under  this  circumstance,  any  spacial  variation  of  the  thickness  could 
be  determined  by  either  an  amplitude  or  phase  measurement  of  the  temperature  at 
the  surface  of  the  film. 

2.3.2  Thermal  Frequency  Response 

The  experimental  work  presented  in  chapters  4  and  5  has  been  done  with  optical 
power  modulated  at  frequencies  near  1  GHz.  This  raises  the  question  of  whether  the 
material  can  thermally  respond  at  such  a  high  frequency.  In  this  section,  the  results 
of  an  analysis  of  thermal  frequency  response  will  be  presented.  The  essence  of  this 
analysis  is  that  when  the  dimensions  of  the  thermal  source  are  smaller  than  or 
equal  to  the  thermal  wavelength,  the  thermal  frequence  response  will  not  limit  the 
sensitivity  of  the  photothermal  system.  Intuitively  this  result  can  be  understood 
by  considering  the  fact  that  the  heat  flux  away  from  the  thermal  source  region 
is  proportional  to  the  temperature  gradient.  When  the  source  is  very  small,  the 
gradient  becomes  very  large,  and  the  flux  can  be  enormous.  The  large  heat  flux 
makes  possible  a  rapid  change  in  the  temperature.  This  lack  of  thermal  inertia  is 
equivalent  to  a  frequency  response  which  extends  out  to  very  high  frequencies. 

The  thermal  frequency  response  of  a  system  is  important  because  the  sensitivity 
of  a  photothermal  technique  is  generally  limited  only  by  the  maximum  allowable 
optical  power  that  can  be  deposited  non-destructively  in  the  sample.  Since  negative 
heat  cannot  be  supplied,  there  is  always  as  much  dc  or  average  heating  as  there 
is  ac  heating.  In  general,  the  ac  response  is  less  than  the  dc  response.  Since  the 
photothermal  signal  is  proportional  to  the  ac  temperature  response,  one  would  like 
to  have  an  ac  response  that  is  near  that  of  the  dc  response  so  that  the  upper  bound 
on  average  heating  does  not  limit  the  size  of  the  ac  photothermal  signal. 


To  examine  the  thermal  frequency  response,  the  program  listed  in  the  Appendix 
was  written.  The  program  calculates  the  ac  and  dc  temperature  rise  at  the  surface 
of  a  semi-infinite  solid  when  heated  by  a  modulated  Gaussian  beam  of  light. 
The  variables  include  the  beam  width,  optical  absorption  depth,  and  modulation 
frequency,  as  well  as  the  thermal  conductivity,  density,  and  specific  heat  of  the 
solid. 

In  Figure  2-3,  the  ratio  of  the  ac  to  dc  response  is  plotted  as  a  function  of  two 
parameters,  the  ratio  of  the  thermal  wavelength  to  optical  absorption  depth  and  the 
ratio  of  the  thermal  wavelength  to  the  optical  spot  radius.  The  plot  indicates  that 
the  thermal  source  should  have  dimensions  smaller  than  or  equal  to  the  thermal 
wavelength  if  the  ac  response  is  to  be  near  the  dc  response. 

Figure  2-4  shows  the  ac/dc  response  ratio  for  1  GHz  modulation  frequency 
in  silicon  as  a  function  of  optical  spot  size.  As  can  be  seen  from  the  plot,  the  ac 
response  can  be  near  the  dc  response  at.  a  frequency  of  1  GHz,  if  the  optical  beam  is 
tightly  focused  to  a  submicron  spot.  Practically,  this  means  that  if  the  silicon  can 
tolerate  an  average  temperature  of  20  degrees  above  the  ambient,  then  a  5  degree 
harmonic  response  at  a  frequency  of  1  GHz  can  be  obtained  if  a  submicron  optical 
spot  is  used.  In  the  next  section,  it  will  become  clear  that  having  a  thermal  source 
which  is  smaller  than  or  equal  to  the  thermal  wavelength  is  important  for  another 
reason. 

2.3.3  Dependence  upon  Thermal  Conductivity 

Using  the  program  found  in  the  Appendix,  the  dependence  of  the  harmonic 
thermal  response  on  thermal  conductivity  was  calculated  for  several  optical  spot 
diameters.  The  calculation  was  done  using  the  material  values  of  silicon  at  a 
modulation  frequency  of  l  GHz.  The  value  for  the  thermal  conductivity  of  silicon 
was  varied  between  25  and  500  \V/m°K.  The  ac  thermal  response  for  each  spot 
size  was  calculated  and  normalized  to  show  the  increasing  sensitivity  to  changes  in 
thermal  conductivity  with  decreasing  spot  size.  Since  many  of  the  photothermal 
techniques  can  measure  both  the  amplitude  and  phase  of  the  photothermal  response, 
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Figure  2-3.  The  nc  to  do  temperature  ratio  at  the  surface  of  a  solid  created 
by  the  absorption  of  a  harmonically  modulated  Gaussian  beam. 
The  ratio  is  plotted  against  two  parameters,  the  ratio  of  ther¬ 
mal  wavelength  to  opioal  spot  radius  and  the  ratio  of  thermal 
wavelength  to  optical  absorption  depth. 


these  were  plotted  independently.  The  dependence  of  the  amplitude  on  thermal 
conductivity  is  shown  in  Figure  2-5.  The  dependence  of  the  phase  is  shown  in 
Figure  2-6. 

The  plots  clearly  show  that  there  is  another  reason  for  focusing  the  optical 
spot  as  tightly  as  possible.  When  the  thermal  source  becomes  large  compared  to 
the  thermal  wavelength,  the  dependence  of  the  harmonic  response  on  the  thermal 
conductivity  goes  to  zero.  When  the  thermal  source  is  smaller  than  or  of  the 
same  dimension  as  the  thermal  wavelength,  the  photothermal  signal  can  be  used  to 
directly  measure  the  thermal  conductivity. 

2.4  Applications  of  High  Frequency  Photothermal  Characterization 

In  the  previous  sections,  it  has  been  shown  that  the  photothermal  response 
can  be  used  to  measure  the  geometrical  and  thermal  properties  of  materials.  In 
particular,  it  was  shown  that  the  photothermal  response  depends  upon  the  thermal 
conductivity  of  the  material  under  evaluation.  In  this  section,  the  areas  where  such 
information  could  be  useful  will  be  examined. 

Because  the  thermal  conductivity  is  a  strong  function  of  the  crystallinity  of  a 
semiconducting  solid,  photothermal  characterization  could  be  used  to  evaluate  the 
crystalline  defect  density  in  these  materials.  For  example,  the  thermal  conductivity 
of  single  crystal  silicon  at  room  temperature  is  near  150  W/m  °  K  [2.5].  Amorphous 
silicon  has  a  thermal  conductivity  of  approximately  2  W/m  0  K  (2.6).  This  is  a 
change  of  almost  two  orders  of  magnitude.  When  a  crystalline  material  is  implanted 
with  high  energy  dopants,  the  lattice  can  be  driven  toward  an  amorphous  state. 
Under  such  circumstances,  high  frequency  photothermal  analysis  could  provide  a 
measure  of  the  resultant  crystallinity  with  high  spacial  resolution. 

When  dopants  are  diffused  into  a  crystalline  semiconductor,  the  thermal 
conductivity  is  changed.  Photothermal  measurement  could  also  provide  a  measure 
of  the  dopant  distribution  and  density.  Holland  [2.7]  has  shown  that  this 
change  is  small  at  room  temperature  when  the  dopant  density  is  less  than 
1018/cm3.  Therefore,  at  room  temperature,  only  a  very  large  dopant  density  can 
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5.  This  plot  shows  the  dependence  of  the  ac  phototherma!  response 
on  the  thermal  conductivity  fur  several  spot  diameters  using  the 
material  values  of  silicon.  The  plot  indicates  that  the  dependence 
on  the  thermal  conductivity  increases  with  decreasing  spot  size. 
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Figure  2-8.  This  plot  shows  the  dependence  of  the  phot  ot  hernial  phase  upon 
the  thermal  conductivity  for  various  spot  diameters  using  the 
material  values  of  silicon.  The  variation  of  the  phase  with  changing 
thermal  conductivity  is  greatest  for  small  spot  diameters. 


be  detected.  However,  at  temperatures  below  100  0  K,  the  thermal  properties 
depend  more  strongly  on  dopant  or  defect  densities  [2.8J.  At  these  temperatures,  the 
photothermal  response  could  provide  a  very  sensitive  measure  of  dopant  densities 
much  less  than  1018/cm3  with  high  spacial  resolution. 

In  anisotropic  crystalline  solids,  the  magnitude  of  the  thermal  conductivity 
depends  upon  the  direction  of  the  diffusion  relative  to  the  primary  crystalline  axes. 
This  fact  manifests  itself  when  a  photothermal  image  of  a  polycrysblline  material 
with  anisotropic  crystallites  is  obtained.  When  the  crystallites  are  oriented  at 
different  angles  with  respect  to  the  focused  beam,  the  photothermal  response  varies. 
This  response  can  be  used  to  identify  and  map  the  various  crystalline  orientations 
in  a  polycrystalline  material. 

As  was  shown  in  section  2.3,  another  application  for  high  frequency  thermal 
waves  is  the  evaluation  of  subsurface  interfaces.  Since  thermal  waves  can  propagate 
through  optically  opaque  materials,  they  provide  a  unique  tool  for  subsurface 
analysis.  In  a  layered  system,  a  variation  of  the  thermal  properties  of  the  materials 
at  an  interface  will  modify  the  reflection  coefficient  seen  by  the  incident  thermal 
waves.  They  can  therefore  be  used  to  probe  the  material  properties  at  the  interface. 
An  example  of  one  interesting  case  for  study  is  the  interaction  that  takes  place 
when  gold  is  deposited  upon  gallium  arsenide  [2.9].  At  the  interface  between  these 
materials,  a  time  dependent  chemistry  takes  place  which  can  modify  the  properties 
of  the  system.  Thermal  waves  analysis  could  provide  useful  information  about  the 
spacial  and  temporal  development  of  such  chemistry. 

Also  shown  in  section  2.3  is  the  fact  that  the  thickness  and  thermal  properties  of 
a  thin  film  can  be  determined  by  photothermal  measurement.  At  high  frequencies, 
the  thermal  wavelength  in  insulators  can  be  less  than  200  A.  This  would  make 
possible  the  measurement  of  thin  films  with  thicknesses  less  than  100  A  with  micron 
lateral  resolution  and  good  precision. 

Finally,  the  photothermal  technique  can  be  used  to  probe  the  adhesion 
properties  of  thin  films  or  layered  systems.  The  advantage  of  using  a  high  frequency 
photothermal  technique  over  other  methods  is  that  the  thermal  wavelength  can  be 
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short  compared  to  optical  or  acoustic  wavelengths.  When  a  separation  between  two 
media  exists,  the  reflectivity  due  to  the  separation  will  be  proportional  the  spacing 
to  wavelength  ratio.  With  the  short  thermal  wavelengths,  very  small  spacings  can 
be  detected. 
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Chapter  3. 


Thermoacoustic  Generation 


3.1  Introduction 

The  photothermal  response  described  in  the  previous  chapter  is  difficult 
to  directly  detect  due  to  the  highly  damped  character  of  the  thermal  waves. 
Fortunately,  this  harmonic  thermal  response  couples  to  an  acoustic  wave  through 
the  periodic  expansion  which  also  takes  place  in  a  medium  when  its  temperature  is 
locally  modulated.  The  thermoacoustically  generated  sound  is  relatively  undamped 
and  may  propagate  away  from  the  source  region,  carrying  with  it  information  about 
the  local  thermal  and  acoustic  properties  of  the  solid.  In  order  to  extract  this 
information  from  the  acoustic  power  radiated  from  the  source  region,  a  thorough 
understanding  of  thermoacoustic  coupling  is  required.  In  this  chapter,  the  three 
dimensional  thermoacoustic  field  equation  will  be  derived  and  solutions  for  one  and 
three  dimensional  generation  will  be  presented. 

In  the  first  section,  the  thermoacoustic  field  equation  will  be  derived  for  an 
infinite  isotropic  solid  with  a  spacially  non-uniform  temperature  distribution.  This 
equation  will  then  be  used  to  derive  the  one  dimensional  thermoacoustic  wave 
equation,  and  it  will  be  shown  that  an  error  exists  in  the  traditional  one  dimensional 
thermoacoustic  theory  in  solids  as  presented  by  White  [3.1].  The  new  theory 
predicts  that  the  thermoacoustic  generation  efficiency  is  nine  times  that  predicted 
by  the  previous  theory.  Using  a  Green’s  function  approach,  the  one  dimensional 
problem  will  be  solved  for  arbitrary  acoustic  boundary  conditions.  This  approach 
will  make  the  understanding  of  the  thermoacoustic  coupling  efficiency  at  an  acoustic 
boundary  very  straight  forward. 

Finally,  the  three  dimensional  thermoacoustic  field  equation  will  be  used  to 
obtain  a  solution  for  a  spherically  symmetric  thermal  distribution  in  an  infinite 
isotropic  solid.  The  implications  of  this  theory  will  be  discussed,  and  a  comparison 
made  with  the  predictions  of  one  dimensional  theory. 


3.2  Three  Dimensional  Thermoacoustic  Field  Equation 


The  stress-strain  relationship  for  an  isotropic  solid  with  uniform  temperature 


is  given  by 


l  =  M  & 


where  X.  is  the  acoustic  stress  tensor,  [c]  is  the  stiffness  tensor,  and  S.  is  the  strain 
tensor.  The  presence  of  non-uniform  temperature  produces  an  additional  strain 
which  does  not  contribute  to  the  stress  (3.2].  If  the  following  definitions  are  made, 

Sff  =  strain  due  to  non-uniform  temperature 
Sa  =  strain  due  to  acoustic  forces  (no  temperature) 

£_  =  Si  +  =  total  strain 

then  the  stress-strain  relation  for  an  isotropic  solid  with  a  non-  uniform  temperature 
distribution  is  given  by 

X-(45l  =  M(£-&)  =  M£-M&-  (3.3) 

The  strain  due  to  the  non-uniform  temperature  is  purely  dilational.  In  contracted 
notation  [3.3],  it  is  given  by 

1 

Si  =  QsO  *  ,  (3.3) 


where  o4  is  the  linear  coefficient  of  expansion  and  0  is  the  temperature  difTerenci 
between  actual  temperature  and  some  reference  temperature  of  the  body.  Takin< 
the  product  of  the  isotropic  stiffness  tensor  and  the  thermal  strain  tensor,  we  obtaii 


[c]S&  =  (<*n  +2ci3)ar,0  !  =3 Baa0  \  , 


where  D  is  the  bulk  modulus.  Under  these  conditions,  the  stress-strain  relation  is 
given  by 


Z  =  [c\S.-ZBat0\ 


(l\ 

1 

1 

0 

0 

Vo 


(3.5) 


This  relation  is  called  the  Duhamel-Neumman  law  [3.4].  This  expression  indicates 
that  the  temperature  distribution  cannot  be  looked  at  as  a  simple  generator  of 
strain.  In  fact,  the  stress  and  strain  variables  are  symmetrically  related  to  the 
non-uniform  temperature.  The  relative  amount  of  strain  or  stress  generated  by 
a  thermal  source  is  determined  by  boundary  and  symmetry  conditions,  as  will  be 
discussed  later. 

To  obtain  the  acoustic  field  equation,  the  divergence  of  the  stress  must  be 
equated  to  the  density-acceleration  product  at  each  point. 


V  1  = 


(3.6) 


After  quite  some  algebra,  the  result  is 


V  •  I  =  (X  +  2 /i)V(  V  •  3)  -  /i(V  X  V  X  3)  -  3 BatV0.  (3.7) 

In  an  infinite  medium,  the  thermoacoustic  source  couples  only  to  displacements 
which  are  dilational.  Therefore  any  displacement  which  is  rotational  need  not  be 
considered.  The  curl  of  the  displacement  can  therefore  be  set  equal  to  zero.  This  also 
means  that  the  gradient  of  the  divergence  of  the  displacement  is  just  the  Laplacian 
of  the  displacement.  Under  these  conditions,  the  vector  displacement  equation  with 
a  thermal  source  term  becomes  (for  dilational  displacements) 

(\+2^V2u=p—+3DaS0.  (3.8) 

This  is  the  thermoacoustic  field  equation  for  an  infinite  isotropic  solid  with  a  non- 
uniform  temperature  distribution. 


3.3  One  Dimensional  Thermoacoustic  Generation 


In  this  section,  the  one  dimensional  thermoacoustic  wave  equation  will  be 
obtained  by  substituting  into  equation  (3.8)  a  harmonic  thermal  source  which 
varies  in  only  one  dimension  of  space.  The  results  obtained  differ  from  those 
presented  by  White  (3.1)  in  his  classic  paper  on  thermoacoustic  generation.  The 
source  of  the  discrepancy  will  be  examined.  After  establishing  the  one  dimensional 
thermoacoustic  equation,  a  Green’s  function  approach  will  be  used  to  find  the 
solution.  This  approach  will  clarify  how  the  thermal  distribution  couples  to  an 
acoustic  wave.  This  analysis  will  also  show  that  the  generated  far  field  strain  in  the 
solid  is  proportional  to  the  Fourier  transform  of  the  thermal  distribution  evaluated 
at  the  value  of  the  acoustic  k-vector.  This  analysis  is  applicable  to  an  arbitrary  one 
dimensional  acoustic  boundary  condition.  The  frequency  dependence  of  the  one 
dimensional  thermoacoustic  conversion  efficiency  will  also  be  discussed. 

3.3.1  Thermoacoustic  Wave  Equation 


The  one  dimensional  thermoacoustic  field  equation  can  be  obtained  by  sub¬ 
stituting  a  one  dimensional  thermal  distribution  into  the  three  dimensional  ther¬ 
moacoustic  equation. 

0(x,  y,  z)  =  0(z)  (3.9) 


If  the  gradient  of  this  temperature  distribution  is  substituted  into  equation  (3.8), 
three  scalar  equations  are  obtained.  The  two  equations  in  the  x  and  y  coordinates 
do  not  have  source  terms  and  therefore  can  be  neglected.  The  equation  in  the  z 
coordinate  is  given  by 


(X  +  2/i) 


d2uz 

dz2 


d2u,  or>  dO 
p-  =  3 Da„  —  . 


dt 2 


dz 


(3.10) 


Under  temporally  harmonic  conditions,  the  equation  becomes 


This  is  the  one  dimensional  thermoacoustic  wave  equation. 

The  thermoacoustic  wave  equation  obtained  by  White  is  identical  to  equation 
(3.11)  above,  except  that  the  thermal  source  term  is  a  factor  of  three  smaller.  The 
discrepancy  between  this  result  and  that  of  White  can  be  understood  if  a  comparison 
is  made  between  the  Duhamel-Neumman  law  and  the  relation  which  was  used  to 
derive  White's  result.  The  Duhamel-Neumman  law  provides  a  relation  between  the 
stress  T3,  the  strain  components  Si,  S2,  and  S3,  and  the  temperature  0. 

T3  =  cnSt  +  +  £2)  —  3  Dae0  (3.12) 

Under  one  dimensional  symmetry  in  an  infinite  medium,  the  strain  in  the  lateral 
dimensions  (Si  and  So)  must  be  zero.  This  is  true  because  there  are  no  boundaries  or 
temperature  gradients  to  break  the  lateral  symmetry.  Thus,  the  individual  elements 
of  the  solid  see  uniform  compressional  forces  from  the  later;- i  directions,  and  the 
net  lateral  displacement  is  zero. 

S,  =S2  =  0  (3.13) 

Substituting  this  condition  into  equation  (3.12)  provides  the  following  relation 
between  the  stress,  strain  and  temperature  under  one  dimensional  symmetry, 

r3  =  (X  +  2//)S3  -Wa„0  (3.14) 

where  (X  +  2/t)  =  Cu.  White  derived  the  thermoacoustic  field  equation  using  the 
following  relation. 

T3=(\  +  2fi)S3-Bae0  (3.15) 

Comparing  these  two  equations,  it  is  apparent  that  there  is  a  factor  of  three  missing 
from  the  relation  used  by  White.  A  look  at  the  reference  [3.2]  sighted  by  White  as 
the  source  of  equation  (3.15)  makes  it  clear  that  an  error  was  made  in  transferring 
this  relation  from  the  reference.  Analysis  done  by  Bray  [3.5]  and  Rosencwaig 
[3.6]  also  contain  the  same  error.  Opsal  [3.7]  uses  the  correct  one  dimensional 
thermoacoustic  field  equation  (3.11)  without  any  comment  or  derivation. 

Since  equation  (3.11)  is  a  field  equation,  a  variation  of  the  magnitude  of 
the  source  term  by  a  factor  of  three  means  an  change  in  the  thermoacoustically 
generated  power  by  a  factor  of  nine.  The  thermoacoustic  conversion  efficiency  is 
therefore  is  almost  an  order  of  magnitude  larger  than  that  predicted  by  the  previous 
theory. 


3.3.2  Solution  of  the  Thermoacoustic  Wave  Equation 


In  order  to  obtain  a  solution  to  the  one  dimensional  field  equation,  it  is  most 
convenient  to  transform  the  equation  from  the  displacement  variable  to  acoustic 
stress.  The  reason  for  transforming  is  the  following.  The  total  strain  in  a  body 
created  by  a  non-uniform  temperature  is  the  sum  of  the  thermally  generated  strain 
and  the  acoustic  strain.  The  acoustic  stress,  however,  is  determined  only  by  the 
acoustic  strain  [3.2].  Therefore,  obtaining  an  equation  in  the  acoustic  stress  makes 
possible  the  separation  of  the  thermal  and  acoustic  elements  of  the  problem.  This 
is  particularly  important  when  using  a  Green's  function  approach.  The  acoustic 
Green's  function  represents  the  propagating  element  of  a  given  acoustic  variable.  If 
the  equation  is  in  the  displacement  or  strain  variable,  these  quantities  include  both 
a  thermal  and  acoustic  contribution.  The  thermal  element  does  not  propagate. 
Therefore,  the  acoustic  Green's  function  cannot  be  used  to  describe  that  portion  of 
the  variable.  If  an  equation  is  obtained  in  acoustic  stress,  the  thermal  and  acoustic 
contributions  are  separated  and  application  of  the  Green’s  function  approach  is  very 
straight  forward.  Another  reason  for  transforming  the  equation  to  the  stress  variable 
is  that  acoustic  boundaries  can  be  easily  handled,  since  the  thermal  distribution 
does  not  modify  the  acoustic  boundary  condition  for  stress. 

The  one  dimensional  thermoacoustic  wave  equation  (3.11)  will  now  be  trans¬ 
formed  to  an  equation  in  the  stress  variable,  using  equation  (3.14)  derived  in  the 
previous  section. 

T3  =  (\+2fi)S3-3Baa0  (3.16) 


In  another  form,  this  relation  becomes 


T3  |  3 Dase 
^ii  Cn 


(3.17) 


Since  the  thermoacoustic  wave  equation  derived  earlier  was  in  the  variable  u:,  the 
total  displacement  in  the  z  direction,  we  can  take  the  derivative  of  that  equation 
with  respect  to  z.  and  obtain  an  equation  in  the  variable  S3. 


d2S3  ,2  c  _3 Bafd~0 

o  i  ^(Tc3  a  o 


(3.18) 
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Substituting  the  relation  found  in  equation  (3.17)  into  this  equation,  the  ther 
moacoustie  equation  in  the  stress  variable  is  obtained. 

d2  (T3  .  3 Baae\  .  t9(  7 3  .  3 Baa0\  3Basd29 


There  are  two  identical  terms  on  each  side  of  the  equation  which  cancel,  leaving 


the  following. 


+  k2T3  =  -3k2Ba,0 


(3.20) 


This  is  the  one  dimensional  thermoaeoustic  wave  equation  in  the  stress  variable. 
The  Green’s  function  for  the  equation  above 


+  —  <5(-^  —  ~o)> 


(3.21) 


is  given  bv 


G(c,r0)  =  -!-^,-'-Jo1. 


(3.22) 


The  solution  to  the  wave  equation  (3.20)  is  given  by  the  integral  of  the  product  of 
Green’s  function  G  and  the  source  function  —3k2Baa0  over  the  source  coordinate 

°  f00  1 

T3(z)  =  /  -3fc*£?o,0(ro)  — e>fc-  ^~Zolc/z0  (3.23) 

J  —  co  ~J*a 

When  z  is  chosen  to  be  far  from  the  thermal  source,  the  equation  becomes 

T3(:)  =  — ■  — [  0(zo)e}'k‘:°dzo  (3.21) 

This  equation  indicates  that  the  far  field  stress  amplitude  is  proportional  to  the 
Fourier  transform  of  the  temperature  distribution,  evaluated  at  the  acoustic  k- 


vcctor 


7i(--|  = 

2  J 

/CO 

■CO 


(3-25) 


(3.20) 


This  is  the  far  field  stress  generated  by  a  one  dimensional  harmonic  thermal 
distribution  in  an  infinite  isotropic  elastic  medium. 


3.3.3  Thermoacoustic  Generation  near  an  Acoustic  Boundary 


To  include  the  effects  of  an  acoustic  boundary  near  the  thermal  source 
distribution,  the  method  of  images  may  be  used.  Adding  a  mirror  image  of  the 
thermoacoustic  source  across  the  acoustic  boundary,  with  the  magnitude  and  phase 
of  the  acoustic  reflection  coefficient  at  the  boundary,  provides  an  effective  thermal 
source  without  the  boundary  which  is  equivalent  to  the  original  distribution  with 
the  acoustic  boundary.  This  is  shown  in  Figure  3-1. 


Otfj  =  0(z)  +  r0{-z) 


(3.27) 


The  acoustic  reflection  coefficient  r  is  determined  by  the  acoustic  impedance 
mismatch  at  the  acoustic  boundary. 


Zi  —  Z2 
Zi  +  Z2 


(3.28) 


Z;  is  the  acoustic  impedence  of  medium  i. 

If  this  result  is  substituted  into  equation  (3.25),  the  far  field  stress  is  given  by 

T3U)  =  ~S^A-e-^:r{(0(zo)  +  rO(-z 0))}|fca,  (3.25) 

It  is  a  property  of  Fourier  transforms  that  if 


?{0(z)}  =  ©(At), 


(3.30) 


7{0[~z)}  =  ©(-A-) 


(3.31) 


Therefore,  the  far  field  stress  generated  by  a  thermal  source  0(~o)  near  an  acoustic 
boundary  with  reflection  coefficient  r  is  given  by 


l-i  =  + 


(3.32) 


where 


©(A„)=7{0(c)}L  . 


(3.33) 
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T!ie  formalism  above  makes  possible  the  simple  calculation  of  all  the  inter¬ 
mediate  cases  between  the  constrained  and  free  surface  boundary  conditions.  To  il¬ 
lustrate  the  effect  of  the  acoustic  boundary  conditions  and  frequency  of  modulation 
on  the  thermoacoustic  conversion  efficiency,  the  following  case  will  be  examined. 
Two  semi-infinite  solids  of  different  acoustic  impedance  are  bonded  together.  One 
of  the  solids  is  transparent  and  thermally  non-conducting.  Harmonically  modulated 
light  is  transmitted  through  the  transparent  solid  to  the  interface  where  it  is  ab¬ 
sorbed  in  an  infinitesimally  thin  boundary  layer  of  the  absorbing  solid.  The  thermal 
distribution  created  by  the  harmonic  heating  couples  to  an  acoustic  wave  which 
propagates  away  from  the  interface.  The  photoacoustically  generated  power  is 
measured  in  the  optically  absorbing  solid  far  from  the  thermal  source.  While  keep¬ 
ing  all  other  parameters  constant,  the  conversion  efficiency  is  measured  as  a  func¬ 
tion  of  the  modulation  frequency  and  the  acoustic  impedance  mismatch.  Figure  3-2 
contains  a  family  of  curves  of  the  conversion  efficiency  versus  modulation  frequency 
for  several  values  of  the  acoustic  reflection  coefficient.  The  reflection  coefficient  is 
related  to  the  impedance  mismatch  at  the  acoustic  boundary  by  equation  (3.28). 

There  are  several  interesting  features  of  these  curves.  At  a  free  surface,  the 
conversion  efficiency  is  proportional  to  frequency.  If  a  photoacoustic  experiment  is 
done  in  air,  the  acoustic  reflection  coefficient  r  is  -.99998.  Under  this  condition, 
there  is  strong  reason  for  working  at  high  frequencies.  Electron-acoustic  imaging 
is  done  under  vacuum  and  again  the  far  field  thermoacoustic  conversion  efficiency 
should  increase  linearly  with  frequency. 

These  results  are  based  upon  a  one  dimensional  model.  In  the  final  section  of 
this  chapter,  the  three  dimension  theory  will  be  presented  to  account  for  the  highly 
localized  nature  of  the  thermal  distribution  which  is  most  common  in  practice. 

3.4  Thermoacoustic  Generation  under  Spherical  Symmetry- 

In  this  section,  the  solution  to  the  thermoacoustic  field  equation  for  a  point 
source  of  harmonically  varying  power  will  be  presented.  The  solution  of  three 
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Figure  3-2.  The  relative  photoacoustic  conversion  elliciency  as  a  function  of 
modulation  frequency  and  acoustic  boundary  conditions  for  one 
dimensional  generation  of  sound. 


dimensional  problems  in  solids  in  general  is  rather  intractable.  This  case  was  chosen 
because  it  could  be  solved  in  a  straight  forward  manner.  The  results  are  general, 
however,  in  the  sense  that  when  the  thermal  source  is  small  compared  to  an  acoustic 
wavelength,  the  actual  shape  of  the  source  becomes  unimportant.  The  results  show 
that  the  conversion  efficiency  increases  strongly  with  frequency. 

3.4.1  Thermoacoustic  Equation  in  Spherical  Coordinates 


The  thermoacoustic  field  equation  for  an  infinite  isotropic  solid  with  a  non- 
uniform  temperature  distribution  was  derived  in  section  3.2.  The  result  is  shown 
below. 

(X  +  2  /<)V23  =  p—  +  3J3aaV0.  (3.34) 

To  obtain  the  solution  to  this  equation  for  a  temperature  source  which  is  spherically 
symmetric,  a  transformation  to  spherical  coordinates  is  necessary.  The  symmetry 
requires  that  all  displacements  must  be  radial,  and  all  derivatives  transverse  to  the 
radius  vector  must  vanish.  In  spherical  coordinates,  equation  (3.34)  becomes  the 
scalar  expression  shown  below. 


„  .  d  (  1  d  2 


A  d 2 


+  3Ba,  ~ 
or 


(3.35) 


Integrating  with  respect  to  r  and  differentiating  with  respect  to  time,  the  equation 
becomes 

(X  +  2/i)(~(r2,v))  =  p^  J  Vrdr  +  3Daa~.  (3.3G) 

If  a  velocity  potential  is  defined  such  that 
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or  equivalently, 


equation  (3.3G)  becomes 
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The  constant  is  indistinguishable  from  a  uniform  temperature  change  of  the  entire 
volume,  which  will  not  effect  a  harmonic  excitation  and  solution.  Therefore,  the 
constant  can  be  set  to  zero,  leaving 
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where  Va 2  =  (X  +  2 ft)/p  and  vr  =  d$/dr.  This  is  the  thermoacoustic  field  equation 
for  a  spherically  symmetric  thermal  source  in  an  infinite  isotropic  elastic  medium. 
In  the  next  section,  the  temperature  distribution  generated  by  a  harmonic  point 
source  of  heat  will  be  calculated  for  substitution  into  this  equation. 

3.4.2  Temperature  Distribution  for  a  Periodic  Point  Source 

The  solution  for  a  harmonically  varying  point  source  of  heat  in  an  infinite 
isotropic  solid  has  been  derived  [3.8].  For  a  harmonic  source  of  power  released  from 
a  point  at  the  origin  of  an  infinite  isotropic  medium, 
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the  steady  state  temperature  distribution  is 
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3.4.3  Spherical  Solution 

Substitution  of  equation  (3.42)  into  the  acoustic  field  equation  (3.40)  results  in 
the  following  expression, 
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The  general  solution  to  this  equation  is  given  by 


To  determine  the  coefficient  a  ,  a  boundary  condition  must  be  used.  By  symmetry, 
the  radial  velocity  must  be  zero  at  the  origin.  This  means  that  the  gradient  of  the 
potential  function  must  be  zero  at  the  origin.  Under  these  conditions,  the  total 


solution  is 


/'e~jk*r  e~k,r\ 

*  =  ?Te(— — )  ,3'JG| 

To  examine  the  efficiency  of  photoacoustic  conversion,  a  calculation  of  the  far 
field  acoustic  power  radiated  from  the  source  can  be  made.  Since  the  total  far  field 
power  radiated  is  independent  of  the  distance  from  the  source,  the  calcuation  can 
be  done  at  a  distance  where  the  plane  wave  impedance  can  be  used.  Under  this 
condition,  the  radial  velocity  is  given  by 
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The  plane  wave  impedance  is 


Zo  =  —  =  -pVa. 

Vr 
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Integrating  the  power  density  over  a  spherical  surface  far  from  the  origin,  the  total 
power  radiated  is 
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In  another  form, 
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This  is  an  expression  for  the  total  far  field  acoustic  power  generated  by  a  harmonic 
point  source  of  heat. 
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3.4.4  Frequency  Dependence  of  Thermoacoustic  Coupling 


We  now  consider  the  frequency  dependence  of  the  photoacoustic  conversion 
efficiency.  In  most  solids,  at  frequencies  below  10  GHz, 
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Under  this  condition,  equation  (3.50)  becomes 
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Note  that  the  radiated  power  is  proportional  to  the  square  of  the  modulation 
frequency.  Thus  the  conversion  efficiency  increases  very  fast  with  modulation 
frequency  for  highly  localized  heat  sources. 

Equation  (3.52)  does  not  include  the  effects  of  an  acoustic  boundary  on  the 
conversion  efficiency.  When  the  source  is  near  an  acoustic  boundary,  the  theory  of 
one  dimensional  generation  can  be  used  to  understand  the  frequency  dependence  of 
the  conversion  efficiency.  The  one  dimensional  theory  showed  that  a  free  surface 
increases  the  frequency  dependence  of  the  conversion  efficiency.  This  is  due  to  the 
fact  that  the  reflections  from  the  free  surface  tend  to  cancel  the  forward  propagating 
waves.  At  low  frequencies,  they  cancel  almost  exactly.  At  higher  frequencies,  the 
cancellation  is  less.  Thus,  the  overall  coupling  increases  with  frequency  in  the 
presence  of  the  free  surface.  This  means  that  the  thermoacoustically  generated 
power  radiated  to  the  far  field  goes  as 
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where  N  >  2  when  the  three  dimensional  thermal  source  is  near  a  free  surface. 

It  is  clear  that  when  the  thermal  source  size  is  small  compared  to  the 
acoustic  wavelength,  the  modulation  frequency  should  be  as  high  as  is  possible 
for  optimal  conversion  efficiency.  One  dimensional  theory  does  not  predict  this  in 
the  constrained  case,  and  predicts  only  an  increase  proportional  to  the  frequency 
in  the  free  surface  case.  This  strong  dependence  of  the  conversion  efficiency  on 
frequency  is  one  reason  for  photoacoustic  imaging  at  high  frequencies.  When 
imaging  with  high  resolution,  the  thermal  spot  size  will  always  be  small  compared  to 
an  acoustic  wavelength.  These  results  are  therefore  very  applicable  for  calculating 
the  conversion  efficiencies. 
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Chapter  4. 


HIGH  RESOLUTION  PHOTOACOUSTIC  IMAGING 


4.1  Introduction 

One  of  the  interesting  characteristics  of  photoacoustic  imaging  is  that  light 
can  be  used  to  periodically  deposit  heat  over  a  region  that  is  small  compared  to 
an  acoustic  wavelength.  Under  this  condition,  the  periodic  thermal  distribution 
generates  a  spherically  diverging  sound  wave.  Efficient  detection  of  the  acoustic 
power  in  this  spherical  wave  is  not  possible  with  phase  sensitive  planar  acoustic 
transducers.  The  reason  for  the  inefficiency  is  the  requirement  that  the  phase  of  the 
wave  vary  no  more  than  one  half  wavelength  over  the  aperture  of  the  transducer. 

If  the  transducer  is  placed  near  the  source,  much  of  the  power  falls  within  the 
transducer  aperture,  but  the  spacial  phase  variation  over  the  aperture  greatly 
reduces  the  efficiency  of  transduction.  If  the  transducer  is  placed  far  from  the 
source  distribution,  the  phase  of  the  incident  acoustic  field  is  uniform  across  the 
aperature,  but  little  of  the  total  power  is  collected.  The  solution  to  the  problem 
of  efficiently  detecting  the  diverging  power  is  to  collimate  it  with  an  acoustic  lens. 
This  collimation  can  improve  the  sensitivity  of  the  detection  system  by  orders  of 
magnitude. 

The  concept  of  using  an  acoustic  lens  to  collimate  and  detect  photoacoustically 
generated  sound  was  first  demonstrated  in  1978  by  Wickramasinghe  et  al.  [4.1]. 
An  acoustic  lens  was  used  to  collect  sound  generated  by  a  highly  focused  optical 
beam.  The  sample  was  scanned  relative  to  the  fixed  optical  and  acoustic  lenses  to 
obtain  a  photoacoustic  image.  The  system  used  for  this  initial  demonstration  was 
not  optimal,  however.  The  laser  used  to  generate  the  sound  was  not  modulated 
efficiently,  generating  sound  over  a  frequency  range  of  more  than  5  GHz.  The 
acoustic  lens  used  was  larger  than  necessary,  requiring  that  the  sound  propagate 


through  a  relatively  long  water  path.  The  acoustic  loss  in  the  water  was  reported 
to  be  32  db.  Decause  the  system  was  pulsed,  the  detection  bandwidth  required 
was  approximately  10  MHz.  .4s  a  result  of  the  inefficiency  of  the  generation  and 
detection  system,  very  high  peak  temperatures  were  required  to  obtain  reasonable 
signals.  The  images  indicated  that  the  required  temperatures  were  near  the  melting 
point  of  the  materials  imaged.  Bray  [-1.2J  continued  this  work  and  indicated 
that  even  after  improvements  to  the  system,  temperatures  of  several  hundred 
degrees  were  required  to  obtain  only  modest  signal  to  noise  ratios.  Since  these 
initial  efforts,  several  improvements  have  been  included  in  a  newly  built  system. 
These  improvements  have  increased  the  sensitivity  of  this  technique  by  orders  of 
magnitude,  requiring  only  small  peak  temperatures  and  optical  peak  powers. 

In  this  chapter,  two  photoacoustic  experiments  will  be  described  which  make 
use  of  an  acoustic  lens  for  the  collimation  and  detection  of  high  frequency 
photoacoustically  generated  sound.  In  the  first  experiment,  an  acoustic  lens  is 
used  to  collect  photoacoustic  sound  generated  in  thin  films  on  optically  transparent 
substrates.  In  the  second  experiment,  the  use  of  gallium  as  a  coupling  liquid  has 
made  possible  the  detection  of  sound  generated  in  optically  opaque  materials.  These 
experiments  will  be  described,  and  images  will  be  presented  demonstrating  the 
sensitivity  and  resolution  of  these  techniques. 

4.2  Experimental  Description 

The  experimental  arrangement  is  shown  in  Figure  4-1.  A  modulated  laser 
beam  is  focused  by  an  optical  microscope  objective  through  the  transparent  sample 
substrate  to  a  small  spot  on  the  film  to  be  imaged.  An  acoustic  lens  is  positioned  to 
be  confocal  with  the  optical  lens.  The  modulated  optical  power  photoacoustically 
generates  sound  in  the  region  of  the  optical  spot.  Because  the  optical  spot  size  is 
generally  smaller  than  an  acoustic  wavelength,  the  sound  diverges  from  the  source 
as  a  spherical  wave.  This  sound  propagates  from  the  heated  surface  to  the  acoustic 
lens  through  a  coupling  liquid,  where  it  is  collimated.  The  collimated  sound  then 
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Figure  4-1.  The  experimental  arrangement  for  photoaeoustie  imaging  of  thin 
films  on  optically  transparent  substrates  using  an  acoustic  lens  to 
collect  the  sound  generated  by  the  modulated  optical  power. 
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propagates  to  a  thin  film  zinc-oxide  piezoelectric  transducer,  where  an  RF  signal 
is  generated.  The  RF  signal  is  sent  to  a  low  noise  preamplifier,  after  which  it  is 
further  amplified,  filtered  and  finally  detected  with  a  crystal  detector. 

The  sample  is  scanned  with  respect  to  the  fixed  acoustic  and  optical  lenses  to 
generate  an  image  which  maps  the  photoacoustic  conversion  efficiency  as  a  function 
of  position  on  the  sample.  The  scanning  is  done  at  30  Hz  in  one  dimension,  and  slow 
scanned  in  the  other  dimension.  The  X  and  Y  position  of  the  sample  is  detected 
using  linear  variable  differential  transformers.  The  two  dimensional  position  and  the 
photoacoustic  signal  are  sent  to  a  scan  converter  which  displays  the  photoacoustic 
signal  as  a  function  of  position  on  a  high  resolution  television  screen.  An  image  can 
be  obtained  in  15  seconds  if  a  30  kHz  detection  bandwidth  is  chosen. 

The  alignment  of  the  acoustic  lens  with  respect  to  the  optical  spot  is  very 
tedious,  requiring  a  two  dimensional  search  for  a  small  signal  which  emanates  from 
the  film  over  an  area  of  approximately  one  square  micron.  To  elimanate  one  of  thu. 
degrees  of  freedom  in  the  search,  a  electromechanical  scanner  was  also  built  to  scan 
the  lens  with  respect  to  a  fixed  sample  and  optical  spot.  By  scanning  the  acoustic 
lens  in  one  dimension  at  30  Hz,  the  photoacoustic  response  over  a  line  is  obtained. 
This  photoacoustic  line  response  can  be  displayed  on  a  CRT  screen  in  real  time. 
By  slowly  translating  the  lens  in  the  direction  perpendicular  to  the  fast  scan,  a  one 
dimensional  search  for  the  photoacoustic  spot  can  be  made. 

Modulation  of  optical  power  at  frequencies  near  1  GHz  is  not  a  trivial  task. 
Fortunately,  external  modulation  was  made  unnecessary  by  the  use  of  a  cavity 
dumped  argon  laser.  The  cavity  dumped  argon  laser  consists  of  a  fused  quartz  Bragg 
cell  placed  in  the  folded  cavity  of  a  Spectra-Physics  Model  166  Argon  ion  laser.  The 
arrangement  is  such  that  light  can  be  diffracted  out  of  the  laser  cavity  by  170  MHz 
sound  which  traverses  the  Bragg  cell.  There  are  two  diffracted  components  which 
leave  the  laser  cavity,  one  which  is  upshifted  by  <170  MHz,  the  other  down  shifted 
by  the  same  frequency.  Because  these  two  components  are  spacially  superimposed, 
they  beat  against  each  other  to  produce  a  power  or  intensity  modulation  at  twice 
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this  frequency,  or  940  MHz.  This  argon  laser  is  capable  of  delivering  20  nanosecond 
optical  pulses  at  a  rate  of  up  to  1  MHz  with  peak  powers  near  50  watts. 

A  second  source  of  modulated  optical  power  used  was  a  gallium  aluminum 
arsenide  diode  laser  [4.3].  It  supplies  up  to  5  mW  of  optical  power  in  a  single 
transverse  and  single  axial  mode.  A  circuit  was  built  to  allow  direct  current 
modulation  of  the  device  at  frequencies  near  1  GHz.  The  modulation  efficiency 
was  measured  on  a  sampling  scope.  See  Figure  4-2.  Taking  into  account  the  roll  off 
of  the  PIN  photodiode  [4.4]  used  to  detect  the  optical  modulation,  the  modulation 
depth  is  estimated  to  be  near  100  percent. 

4.3  Sensitivity  of  the  Acoustic  Lens  as  a  Photoacoustic  Detector 

Since  the  photoacoustic  conversion  efficiencies  are  rather  small,  the  sensitivity 
of  the  acoustic  lens  makes  it  attractive  for  collecting  photoacoustically  generated 


Figure  4-2.  The  modulation  depth  of  the  output  of  an  NEC  NDL3000  diode 
laser  at  a  frequency  of  I  GHz. 


power  Using  the  three  dimensional  theory  presented  in  Chapter  3,  a  calculation  can 
be  made  of  the  expected  signal  to  noise  ratio  under  typical  photoacoustic  imaging 
conditions.  This  calculation  is  shown  in  Table  4-1  below. 


TABLE  4-1 


Photoacoustic  Sensitivty  with  an  Acoustic  Lens 


Optical  power  focused  to  1  /i2  area 

io-3 

Watts 

Photoacoustic  conversion  efficiency 

10~9 

Generated  acoustic  power 

10-12 

Watts 

Acoustic  collection  efficiency 

io-2 

Detected  power 

lo-n 

Watts 

Johnson  noise  power  ( 104  Hz  bandwidth) 

10-10 

Watts 

Signal  to  Noise  ratio 

100 

To  experimentally  verify  the  predicted  sensitivity,  a  current  modulated  diode 
laser  [1.3]  was  focused  to  approximately  one  square  micron  onto  the  surface  of  a 
1000  A  thick  gold  film  pattern  on  a  glass  substrate.  The  object  was  then  scanned 
across  the  pattern,  generating  the  a-trace  seen  in  Figure  4-3.  Approximately  0.5 
m\V  of  the  optical  power  was  absorbed  by  the  gold  film.  The  bandwidth  of  the 
detection  system  was  10  kHz.  The  experimentally  measured  sensitivity  is  near  that 
calculated  in  Table  4-1. 

4.4  Photoacoustic  Imaging  Results 

In  this  section,  the  images  obtained  with  the  photoacoustic  imaging  system 
described  in  the  previous  sections  will  be  presented  and  discussed.  The  images 
presented  include  recrystallized  silicon  and  gold  films  on  transparent  substrates. 


Figure  4-3.  The  signal  to  noise  ratio  obtained  using  less  than  1  mW  of  optically 
absorbed  power  to  generate  sound  in  a  1000  A  thick  gold  film  on  a 
glass  substrate. 

4.4.1  Laser  Recrystallized  Silicon  Film  on  a  Quartz  Substrate 

Some  interest  has  been  shown  in  laser  annealing  and  recrystallization  of  thin 
films  deposited  upon  various  substrates  [4. 5, 4. 6].  Laser  recrystallization  of  solids 
is  inherently  traumatic,  requiring  rapid  temperature  changes  and  large  spacial 
gradients  of  temperature.  Furthermore,  the  melting  and  subsequent  recrystallization 
of  thin  films  takes  place  next  to  a  substrate  whose  lattice  constant,  expansion 
coefficient  and  stiffness  are  generally  different  than  those  of  the  film.  These 
conditions  generally  leave  the  recrystaliized  film  cracked  and  under  stress. 

Since  the  photoacoustic  conversion  efficiency  is  integrally  tied  to  the  mechanical 
and  thermal  properies  of  a  material,  it  was  thought  that  photoacoustic  imaging  of 
these  films  might  provide  some  information  about  the  recrystallization  process.  With 
this  in  mind,  a  recrystaliized  silicon  sample  was  obtained  from  the  research  group 


of  Prof.  J.F.  Gibbons  group  at  Stanford  University  [4.7].  Using  an  cw  argon  ion 
laser  with  some  beam  shaping,  a  0.5  micron  polysilicon  film  on  a  quartz  substrate 
with  500  A  grain  size  was  recrystallized  to  produce  grains  with  dimensions  up  to  100 
microns.  Simultaneously,  large  thermally  induced  cracks  in  the  recrystallized  film 
were  also  produced.  These  can  be  clearly  seen  in  the  acoustic  and  photoacoustic 
images  in  Figure  4-4  on  the  following  page.  Another  region  of  the  recrystallized  film 
is  shown  in  Figure  4-5  with  optical  (DIG)  and  acoustic  comparisons. 

Within  these  obvious  cracks,  faint  evidence  of  the  grain  and  subgrain  bound¬ 
aries  can  be  seen  in  the  acoustic  images.  The  contrast  from  grain  to  grain  is  very 
small,  however.  In  the  photoacoustic  images,  there  is  a  strong  contrast  between  the 
grains.  A  measurement  of  the  photoacoustic  contrast  indicated  that  the  variation 
in  photoacoustic  power  from  grain  center  to  grain  center  was  as  large  as  8  db.  This 
is  shown  by  the  photoacoustic  a-trace  in  Figure  4-6. 

It  is  apparent  from  the  images  in  figures  4-4  and  4-5  that  both  the  optical  (DIC) 
and  the  acoustic  images  contain  much  more  information  about  the  topography 
of  the  sample  than  does  the  photoacoustic  image.  The  strong  dependence  upon 
topography  makes  it  difficult  to  extract  materials  information  from  these  optical 
and  acoustic  images. 

The  large  change  in  the  photoacoustic  conversion  efficiency  from  grain  to  grain 
was  somewhat  puzzling.  Using  a  one  dimensional  numerical  routine  for  modeling 
the  photoacoustic  generation  [4.8],  the  far  field  photoacoustic  power  was  calculated 
while  varying  the  material  parameters  of  the  silicon.  The  results  indicated  there 
was  no  reasonable  variation  of  the  material  properties  that  could  explain  the  large 
contrast  seen  in  these  images.  The  calculations  showed  that  the  photoacoustic 
power  generated  was  almost  independent  of  the  thermal  conductivity,  and  acoustic 
stiffness.  There  was  a  strong  dependence  upon  the  expansion  coefficient  and  specific 
heat,  but  it  was  argued  that  these  material  properties  could  not  vary  by  factors  of 
two  from  grain  to  grain.  The  optical  absorption  was  eliminated  as  a  possible  source 
of  contrast  by  imaging  the  sample  in  optical  transmission  both  under  very  strong 


Figure  4-4.  Acoustic  and  photoacoustic  images  of  an  area  of  laser  recrystalized 
silicon  film  on  a  quartz  substrate. 


Figure  4-5.  A  comparison  of  acoustic,  photoacoustic  and  optical  (DIO)  images 
of  the  laser  recrystalizcd  silicon  sample. 


Figure  4-0.  An  a-trace  illustrating  the  variation  of  the  photoacoustic  signal 
from  grain  to  grain  in  the  recrystallized  silicon  sample. 

and  weak  illumination.  These  optical  images  showed  no  contrast  between  the  grains 
which  had  shown  large  photoacoustic  variations. 

A  suggestion  was  made  (4.9]  that  the  source  of  contrast  could  be  the  orientation 
of  the  grains.  Because  the  generation  efficiency  is  dependent  upon  the  bulk  modulus 
of  the  material,  it  seemed  possible  that  the  anisotropy  of  the  silicon  could  make 
the  conversion  efficiency  dependent  upon  the  crystalline  orientation.  A  calculation 
was  therefore  made  of  the  effective  bulk  modulus  of  silicon  as  its  stiffness  tensor  is 
rotated.  This  modulus  was  then  used  in  the  numerical  calculation  of  the  conversion 
efficiency.  The  result  was  that  the  predicted  far  field  photoacoustic  power  did 
depend  upon  the  orientation  of  the  grains.  However,  the  largest  predicted  variation 
in  the  signal  due  to  orientation  was  1.3  db,  much  smaller  than  the  measured 
variation.  Another  confirmation  that  the  source  of  the  contrast  was  not  grain 
orientation  was  some  diffractometer  analysis.  The  analysis  indicated  that  much  less 


than  one  percent  of  the  grains  had  orientations  different  from  <100>  normal  to 
the  surface.  The  contrast  seen  in  the  photoacoustic  images  was  very  wide  spread, 
involving  greater  than  20  percent  of  the  sample  area.  Therefore,  orientation  is  ruled 
out  as  the  source  of  contrast. 

An  alternative  source  of  contrast  is  the  adhesion  of  the  grains  to  the  substrate. 
If  the  grains  are  not  uniformly  bonded  to  the  substrate,  the  photoacoustic  conversion 
efficiency  will  vary  from  point  to  point.  This  possibility  was  also  numerically 
modeled.  The  results  indicated  that  the  expected  variation  of  the  conversion 
efficiency  from  a  well  adhered  grain  to  a  grain  that  was  completely  free  was  5.5 
db.  This  is  the  most  plausible  explanation  of  the  contrast  seen  in  the  photoacoustic 
images. 

Confirmation  of  this  conclusion  came  in  two  ways.  First,  the  acoustic  reflection 
images  showed  high  reflectivity  in  the  regions  where  the  photoacoustic  conversion 
efficiency  was  low.  This  is  consistent  with  the  adhesion  model  because  when  a  grain 
is  poorly  adhered  to  the  substrate,  the  acoustic  reflectivity  in  that  region  should 
be  unity.  The  acoustic  contrast  between  well  and  poorly  bonded  grains  was  small, 
however.  This  is  due  to  the  fact  that  the  impedance  mismatch  between  w'ater  and 
the  silicon  is  very  large.  Even  when  a  grain  is  well  adhered  to  the  substrate,  there  is 
a  very  large  acoustic  reflection  at  the  silicon/water  interface.  Second,  application  of 
scotch  tape  to  the  sample  surface  did  remove  some  of  the  dark  grains  after  repeated 
attempts.  The  grains  were  difficult  to  remove  how’ever.  This  could  be  attributed  to 
the  very  sharp  boundaries  that  existed  between  grains.  Under  these  conditions,  it 
is  conceivable  that  even  with  no  bonding  to  the  substrate,  the  grains  could  still  be 
difficult  to  extract  from  their  locations.  These  results  led  to  the  conclusion  that  the 
source  of  contrast  in  the  images  is  due  to  the  adhesion  properties  of  the  indiviual 
grains. 

4.4.2  Thin  Gold  Film  on  a  Glass  Substrate 

A  patterned  1000  A  gold  film  on  a  glass  substrate  was  imaged  using  the 
photoacoustic  imaging  system  described  in  section  4.2.  The  results  are  shown  in 


Figure  4-7.  The  optical  image  (c)  indicates  that  there  are  small  babbles  or  defects 
below  the  film  in  certain  areas.  In  the  photoacoustic  image  (a),  it  is  very  apparent 
that  in  these  regions  the  photoacoustic  conversion  efficiency  is  very  much  reduced 
relative  to  the  average  photoacoustic  response.  As  was  shown  in  the  previous  section, 
the  photoacoustic  conversion  efficiency  is  strongly  dependent  upon  the  adhesion  of 
films  to  the  substrate.  This  is  the  source  of  contrast  in  the  photoacoustic  image. 
Only  a  slight  contrast  is  seen  in  the  acoustic  image.  This  is  due  to  the  fact  the 
the  acoustic  lens  is  focused  onto  the  surface  of  the  film.  At  this-  focal  position,  the 
acoustic  sensitivity  to  adhesion  is  rather  poor. 

The  photoacoustic  technique  is  inherently  more  sensitive  to  the  adhesive 
properties  of  thin  films  than  acoustic  reflection  for  the  following  reason.  If  the 
thickness  of  a  film  is  very  small  compared  to  an  acoustic  wavelength,  it  only 
slightly  perturbs  the  acoustic  impedance  and  velocity  of  the  substrate.  Since  the 
mismatch  in  impedance  between  water  and  most  solids  is  very  large,  the  reflection 
from  the  substrate  is  almost  independent  of  the  properties  of  the  film.  This  large 
background  reflectivity  limits  the  sensitivity  of  acoustic  reflection  techniques  when 
examining  thin  films  on  substrates.  In  contrast,  the  photoacoustic  signal  does  not 
contain  a  large  background  signal.  The  boundary  condition  under  which  sound 
is  photoacoustically  generated  can  alTect  the  conversion  efficiency  by  orders  of 
magnitude  as  was  shown  in  section  3.3.3.  This  makes  it  an  attractive  technique 
for  examination  of  thin  films  with  adhesion  problems. 

4.5  Photoacoustic  Imaging  of  a  Free  Surface 

One  limitation  of  the  technique  just  described  is  the  requirement  that  the 
sample  must  have  a  transparent  substrate.  The  reason  for  this  is  that  the  imaging 
is  done  in  transmission.  The  optical  power  is  focused  from  one  side  and  the  acoustic 
lens  collects  the  sound  from  the  other.  Since  the  velocity  and  impedance  mismatch 
between  water  and  most  solids  is  very  large,  the  acoustic  lens  must  be  focused 
directly  onto  the  surface  to  be  imaged  while  the  optical  beam  is  focused  to  the 
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ure  4-7.  A  comparison  of  photoacoustic,  acoustic  and  optical  images  of 
a  gold  film  on  a  glass  substrate. 


There  are  many  important  materials  and  systems  which  do  not  have  optically 
transparent  substrates.  These  materials  cannot  be  evaluated  with  the  previous 
technique.  While  these  material  may  be  optically  opaque,  they  generally  are 
acoustically  transparent  well  beyond  1  GHz  frequencies.  It  is  therefore  conceivable 
to  focus  the  optical  power  onto  the  front  surface  of  these  materials  and  collect  the 
sound  through  the  substrate  or  bulk  of  the  sample.  Figure  4-8  contains  a  comparison 
of  the  two  configurations.  The  upper  figure  is  the  water  coupled  system  described 
in  section  4.2.  In  the  lower  figure,  the  sample  has  been  reversed  to  allow  the  light 
to  be  directly  focused  on  to  the  surface  to  be  imaged.  The  photoacoustic  power  is 
then  collected  through  the  sample  substrate  or  bulk. 

There  are  several  reasons  why  collecting  the  sound  through  the  sample 
substrate  is  difficult.  Since  there  is  a  large  acoustic  impedance  mismatch  between 
solids  and  most  liquids,  the  power  transfer  across  the  interface  can  be  very  poor. 
The  velocity  mismatch  between  liquids  and  solids  causes  considerable  aberration  in 
the  focused  acoustic  beam,  further  reducing  the  collection  efficiency.  The  strong 
refraction  at  the  interface  also  makes  it  difficult  to  focus  the  lens  very  far  into  the 
interior  of  the  solid. 

These  problems  were  addressed  by  Jipson  [4.10]  in  1979.  lie  demonstrated  the 
use  of  an  acoustic  lens  in  imaging  interior  planes  of  a  solid  using  liquid  gallium  as 
a  coupling  liquid.  Liquid  gallium  has  three  important  features.  First,  it’s  acoustic 
impedance  is  very  near  that  of  most  solids.  Therefore  the  power  transfer  across 
the  solid/gallium  interface  is  good.  Second,  the  acoustic  velocity  of  the  gallium  is 
approximately  two  times  larger  than  that  of  water.  This  means  that  the  refraction 
that  takes  place  at  the  interface  between  the  solid  and  the  gallium  is  much  more 
tolerable.  The  third  feature  is  a  very  low  attenuation.  This  is  important  because 
large  acoustic  lenses  are  required  to  image  interior  planes  of  a  solid,  since  the 
backsurface  of  the  object  must  be  brought  into  focus  while  the  lens  remains  some 
finite  distance  away  from  the  front  surface. 

The  work  of  Jipson  indicated  that  collection  of  photoacoustically  generated 
sound  through  the  bulk  of  the  sample  could  be  done.  To  demonstrate  this  possibility, 
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Figure  4-8.  A  comparison  of  the  two  ways  in  which  photoacoustically  generated 
sound  may  be  collected  by  an  acoustic  lens.  In  the  upper  figure, 
the  sound  is  collected  from  the  surface  of  the  solid  with  with  water 
coupling.  In  the  lower  figure,  collection  takes  place  through  the 
hulk  of  the  sample  with  gallium  as  the  acoustic  coupler  between 
the  sample  and  the  lens. 


several  large  acoustic  lenses  were  made  [4.11].  The  radius  chosen  was  250  microns, 
which  allowed  imaging  through  150  microns  of  solid  with  some  reasonable  working 
distance.  The  lenses  were  tuned  to  940  MHz.  The  same  laser,  scanning  system, 
and  RF  detection  system  were  used  as  with  the  case  of  the  water  coupled  system 
described  in  section  4.2. 

A  sample  was  prepared  by  depositing  a  patterned  1000  A  gold  film  upon  a  140 
micron  glass  substrate.  The  sample  was  then  placed  in  the  photoacoustic  imaging 
system  with  the  film  toward  the  optical  lens.  The  film  was  first  imaged  through 
the  substrate  in  pulsed  acoustic  reflection  mode.  This  is  shown  in  the  left  side  of 
Figure  4-9.  The  photoacoustic  image  is  shown  on  the  right. 

The  resolution  of  the  acoustic  image  is  poor  relative  to  that  seen  in  the 
photoacoustic  image.  There  are  two  reasons  for  this  fact.  The  first  is  that  since 
the  acoustic  velocity  of  the  solid  is  high,  the  wavelength  of  the  incident  sound  is 
large  compared  to  the  optical  spot  size  which  essentially  determines  the  resolution 
in  the  photoacoustic  case.  The  second  is  that  the  refraction  that  takes  place  at  the 
solid  gallium  interface  also  spreads  the  acoustic  focus.  These  two  effects  combine 
to  produce  a  poorer  resolution. 

It  is  worth  noting  that  when  imaging  in  acoustic  reflection,  the  acoustic 
reflection  coefficient  as  seen  from  the  substrate  is  near  unity  whether  or  not  the 
patterned  gold  film  is  present  on  the  substrate.  The  acoustic  contrast,  therefore,  is 
largely  due  to  the  scattering  of  the  focused  sound  at  the  boundaries  of  the  gold  film 
rather  than  an  impedance  variation  in  the  sample.  For  this  reason,  the  contrast 
seen  in  the  acoustic  image  of  Figure  4-9  is  largely  found  at  the  edges  of  the  gold 
film  pattern. 

The  photoacoustic  imago  indicates  micron  resolution.  There  is  an  interesting 
interference  elTect  in  the  image.  Apparently,  the  localized  thermal  excitation  couples 
to  a  surface  wave  which  propagates  to  the  film  boundary  where  it  is  partially 
reflected.  It  then  propagates  back  to  the  source  region,  reradiating  some  of  its 
energy  out  into  the  solid.  The  reradiated  energy  interferes  with  the  thormoacoustic 


gure  4-9.  A  comparison  of  acoustic  and  photoacoustic  images  of  a  gold 
film,  obtained  using  gallium  to  collect  the  sound  through  the 
glass  substrate  of  the  sample. 


source  radiation  producing  a  spacially  periodic  variation  of  the  signal  near  the 
boundaries  of  the  patterned  gold. 

The  wetting  of  the  acoustic  lens  with  the  liquid  gallium  was  not  found  to  be 
very  difficult  if  the  lens  and  gallium  were  kept  clean.  The  fact  that  the  radius  of 
curvature  of  the  lens  was  relatively  large  (250  /jm)  helped  minimize  the  wetting 
problems.  Mercury  was  also  used.  The  difficulty  found  in  using  mercury  is  that 
its  surface  tension  is  too  large.  The  only  way  to  use  it  is  to  obtain  a  mercury  ball 
somewhat  smaller  than  the  hemispherical  lens  and  place  it  inside.  However,  with 
the  sample  moving  rapidly  below  the  ball  of  mercury,  it  is  difficult  to  keep  the  ball 
inside  the  acoustic  lens  for  long  periods  of  time.  It  did  seem  possible,  however,  that 
if  a  large  lens  was  designed  to  more  completely  enclose  the  mercury,  the  mercury 
could  be  a  good  coupler  from  the  lens  to  the  sample. 

This  technique  was  also  used  to  image  the  recrystaDi  zed  silicon  sample  discussed 
in  section  4.4.  In  order  to  achieve  proper  wetting  of  the  gallium  on  the  silicon,  it 
is  very  critical  that  the  silicon  surface  be  very  clean.  The  images  obtained  using 
this  technique  also  confirmed  the  non-uniform  adhesion  of  the  grains  in  the  laser 
recrystallized  silicon  film. 
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Chapter  5. 


HIGH  RESOLUTION  PHOTOTHERMAL  LASER  PROBE 


5.1  Introduction 

At  frequencies  near  1  GHz,  the  acoustic  wavelength  in  argon  gas  becomes 
half  that  of  visible  light,  making  possible  collinear  Bragg  scattering  of  laser 
light  by  the  sound.  This  acousto-optic  interaction  has  been  exploited  in  a  new 
phototliermal  laser  probe  to  provide  a  very  high  resolution  technique  for  measuring 
the  photothermal  response  of  materials.  The  probe  has  several  features  which  are 
attractive.  Since  the  operating  frequency  is  high,  the  thermal  wavelength  in  solids 
is  generally  small  compared  to  the  wavelength  of  light.  Under  this  condition,  the 
resolution  is  determined  by  the  optical  spot  diameter  of  the  probe  beam,  which  can 
be  less  than  one  micron.  The  probe  is  an  optical  technique  which  does  not  require 
either  acoustic  transducers  or  contact  with  the  solid  being  evaluated.  Only  small 
optical  input  powers  are  necessary  to  obtain  reasonable  signal  to  noise  ratios  due 
to  its  high  sensitivity. 

The  photothermal  probe  has  been  implemented  in  the  following  arrangement. 
An  intensity  modulated  pump  beam  is  tightly  focused  onto  the  surface  of  a  solid 
in  a  cell  filled  with  argon  gas  at  high  pressure.  See  Figure  5-1.  The  modulated 
pump  beam  generates  a  harmonic  temperature  distribution  which  couples  to  an 
acoustic  wave  in  the  gas  at  the  surface  of  the  solid.  A  second  beam,  the  probe, 
is  collinearly  aligned  with  the  pump  and  also  tightly  focused  to  the  same  location 
on  the  sample.  When  the  modulation  frequency  of  the  pump  beam  is  near  1  GHz, 
the  collinear  Bragg  scattering  condition  is  met  and  a  portion  of  the  probe  beam 
is  back  scattered  or  back  reflected  from  the  sound  propagating  in  the  gas.  The 
light  which  is  scattered  from  the  sound  in  the  gas  is  upshifted  in  frequency,  so  it 
can  be  heterodyned  in  an  optical  detector  with  the  probe  light  which  is  reflected 


Figure  5-!.  The  geometry  of  the  photothermal  laser  probe. 
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from  the  solid  surface.  See  Figure  5-2.  As  will  be  explained  in  later  sections,  the 
heterodyne  signal  is  proportional  to  the  ac  surface  temperature  over  the  region  of 
the  solid  illuminated  by  the  tightly  focused  probe  beam.  A  photothermal  image  is 
obtained  by  scanning  the  object  in  two  dimensions  relative  to  the  pump  and  probe 
beams.  This  provides  a  measure  of  the  thermal  properties  of  the  solid  with  high 
spacial  resolution. 

In  this  chapter,  the  theoretical  basis  for  the  photothermal  sensitivity  of  the 
probe  will  be  presented.  The  sensitivity  will  be  calculated  with  argon  gas  as  the 
coupling  medium.  A  figure  of  merit  will  also  be  established  by  which  various 
coupling  media  may  be  evaluated.  A  comparison  will  then  be  made  between 
three  high  resolution  photothermal  and  photoacoustic  techniques.  Finally,  an 
experimental  description  and  verification  of  the  theory  will  be  presented  including 
high  resolution  photothermal  images  of  boron  implanted  silicon. 


5.2  Derivation  of  the  Bragg  Scattering  Efficiency 

Acoustic  Bragg  scattering  of  light  is  a  non-linear  interaction.  It  can  be 
treated  formally  by  calculating  the  non-linear  polarization  created  when  optical 
and  acoustic  fields  are  superimposed.  This  non-linear  optical  polarization  radiates 
light  at  the  sum  and  difference  frequency  of  the  two  interacting  optical  and  acoustic 
fields.  The  total  power  reradiated  at  the  sum  and  difference  frequency  depends  upon 
the  summation  of  the  fields  over  the  source  region.  In  this  section,  it  will  be  shown 
that  the  Bragg  scattered  electric  field  strength  is  proportional  to  an  integral  of  the 
non-linear  source  polarization  over  the  interaction  region.  The  strength  of  the  non¬ 
linear  polarization  will  then  be  calculated  using  the  photoelastic  relation  for  the  gas. 
Finally,  these  results  will  be  combined  to  provide  a  relationship  between  the  driving 
acoustic  field  in  the  gas  and  the  scattered  electric  field  strength.  From  this  relation, 
an  effective  Bragg  reflectivity  as  seen  by  the  probe  beam  will  be  established. 
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Figure  5-2.  An  illustration  of  the  two  reflected  components  of  the  incident 
probe  beam.  One  component  is  reflected  from  the  surface  of 
the  solid,  while  the  other  component  is  reflected  and  upshifted 
in  frequency  by  the  sound  propagating  in  the  gas. 


5.2.1  Inhomogeneous  Wave  Equation  with  Acousto-optic  Source  Term 


Starting  with  Maxwell’s  equations  and  the  constitutive  relations  for  an  isotropic 
lossless  medium,  one  obtains  the  familiar  one  dimensional  traveling  wave  equation 
with  a  source  polarization  term  (5.1). 

d2£  ,  n2  d2£  d2PNL 

dz 2  +  e2  dt2  ~  /l  dt2  *5 

£  is  the  scattered  electric  field  strength  as  a  function  of  time  and  space,  n  is  the 
index  of  refraction  at  the  frequency  of  interest,  c  is  the  velocity  of  light  in  vacuum, 
/<  is  the  permiability,  and  PNL  is  the  nonlinear  polarization,  also  a  function  of 
time  and  space.  The  scattered  electric  field  component  £  and  nonlinear  source 
polarization  PNL  are  defined  to  be  of  the  following  form. 


£  =  £(z,  l)  =  Re{Eej{ut~kz)}  (5.2) 


PSL  =  P*L(z,t)  =  f?e{P,vV(wt-fcs)} 


(5.3) 


E  and  PNL  are  envelope  functions  of  £  and  P^L.  Both  are  functions  of  time  and 
space.  If  these  definitions  are  substituted  into  the  wave  equation  (5.1),  and  the 
slowly  varying  envelope  approximation  is  employed,  the  result  is 


dE  n  dE  .  1  nNI 
T  +  —57  =  -}-uiiPNL. 
dz  c  dt  2 


(5.4) 


In  steady  state,  dE/dt  =  0,  which  leaves 


dE 

dz 


-i\u,,pNL 


(5.5) 


and 

E  =  i  (5.6) 

This  is  the  solution  for  the  envelope  function  of  the  scattered  electric  field  as  a 
function  of  the  nonlinear  polarization  envelope  function. 


5.2.2  Derivation  of  the  Photoelastic  Relation  for  an  Ideal  Gas 


Before  calculating  the  nonlinear  source  polarization  term,  the  photoelastic 
relation  in  the  gas  must  be  derived.  It  can  be  obtained  using  the  Lorentz-Lorenz 
formula  (5.2) 


3e0(*  -  1) 

°poi  -  7V(c  +  2)  ’ 


(5.7) 


where  o poi  is  the  mean  polarizability  of  the  atoms  or  molecules,  e0  is  the  permittivity 
of  free  space,  e  is  the  relative  dielectric  constant, and  N  is  the  number  of  atoms  or 
molecules  per  unit  volume.  For  a  gas,  where  (c  —  1)  is  much  less  than  1,  it  can  be 
shown  that 


(<-»» 

M\V€o 


(5.8) 


where  MW  is  the  molecular  or  atomic  weight  of  the  medium,  and  p  is  the  density. 
The  derivative  of  equation  (5.8)  provides  a  relation  between  a  change  in  density  of 
the  gas  bp  and  the  corresponding  change  in  the  relative  dielectric  constant  bt.  The 


result  is 


be  =  (e-  l)^  =  (e-l)S, 
Po 


(5.9), 


where 


S  =  £ — £l  (5.10) 

Po 

is  the  acoustic  condensation.  Equation  (5.9)  is  the  photoelastic  relation  for  an  ideal 


gas. 

5.2.3  Calculation  of  the  Acousto-optic  Nonlinear  Polarization 


The  nonlinear  polarization  PNL  can  now  be  calculated  using  the  constitutive 
and  photoelastic  relations  in  the  gas,  and  the  two  incident  optical  and  acoustic 
fields,  So  and  S.  The  constitutive  relation  in  the  gas  is  given  by 


The  change  in  polarization  bP  caused  by  a  change  in  the  dielectric  constant  be, 
given  an  electric  field  £ ,  is 

bP=(0be£.  (5.12) 

The  relation  between  acoustic  condensation  S  and  the  perturbation  of  the  dielectric 
constant  has  already  been  calculated  in  the  previous  section. 

be  —  (e  —  1)5  (5.13) 

Combining  these  results,  the  nonlinear  polarization  can  be  calculated. 

PNL  =  bP  =  e0(e-l)S£.  (5.14) 

The  nonlinear  polarization  is  thus  proportional  to  the  product  of  the  incident 
acoustic  and  optical  fields. 

The  counterpropagating  acoustic  and  optical  driving  fields  are  defined  as 
follows, 

S(z,  t)  =  Re{S(z)ej{w-t~k-z)}  (5.15) 

£0(z,  t )  =  /?e{E0(-y)ej(u'ot+fco2)}  (5.16) 

where  5(c)  and  Eo(z)  are  the  slowly  varying  envelope  functions.  See  Figure  5-3. 
Substitution  of  these  definitions  into  equation  (5.14)  provides  an  expression  for  the 
non-linear  source  polarization. 

P*L  =  if0(f  -  l)^5o(c)Eo(r)eJ'l(u-0+u'')t+^°-fc')l>  +c.c. 

+  50(2)*Eo(c)eJ^a'0-u',^+^'0+fc*^I  +  c.c. 

(5.17) 

Previously,  this  polarization  PNL  has  been  defined  as 

PNL  =  PNL{s,t)  =  Re{PNLe1'{wt-kz)}.  (5.18) 

Equating  the  two  expressions  requires  that  the  envelope  function  PXL  be  given  by 
pnl  =  Ieo(f  _  i)^50(-)zr0(.)eJ'l(“’o+“'.-^)t+(^-fc.+fc)^l 

+  5(2)*£,0(j)eJ'l(tt’°~u'*-u;lt+fA;o+A:*+fclfl 

(5.19) 
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Figure  5-3  The  collinear  Bragg  scattering  geometry. 


5.2.4  Bragg  Scattering  Efficiency 

The  expression  for  the  envelope  of  the  nonlinear  source  polarization  PSL  can 
now  be  used  to  calcualate  the  strength  of  the  scattered  electric  field.  Recall  from 
section  5.2.1  that  the  scattered  electric  field  is  given  by  an  integral  of  the  nonlinear 
source  polarization  envelope  over  the  interaction  region. 

E(=,  t)  =  Jo  (-^r/P  v V,  /))</-'  +  £-(0)  (5.20) 

Substitution  of  the  expression  for  the  nonlinear  polarization  envelope  (5.19)  into 
this  integral  results  in  the  following, 


where  5(c)  has  been  chosen  to  be  real  This  choice  does  not  reduce  the  generality 
of  the  solution. 

The  two  terms  in  the  integrand  will  contribute  little  to  the  integral  if  the 
exponentials  are  oscillating  rapidly  over  the  slowly  varying  source  polarization 
distribution.  This  means  that  for  either  term  to  significantly  contribute,  the  k- 
vector  sum  in  the  exponential  must  be  close  to  zero.  This  is  what  is  termed  as 
k-vector  matching.  See  Figure  5-4. 

This  k-vector  matching  is  not  possible  for  the  second  term  because  the 
magnitude  of  A  would  have  to  be  larger  than  the  magnitude  of  kQ,  but  the  magnitude 
of  u  would  have  to  be  smaller  than  the  magnitude  of  o-’o.  This  cannot  occur  since 
the  k-vector  magnitude  is  essentially  proportional  to  frequency  when  there  is  only 
weak  dispersion.  Thus,  the  second  term  does  not  contribute  significantly  to  the 
integral.  If  the  k-vector  mismatch  if  defined  as 

AA  =  A0  -  A,  -I-  A,  (5.22) 

then  the  expression  for  the  envelope  of  the  scattered  electric  field  is  given  by 

£(c,  /)  =  ~  [  S(z')E0(z')^^k:'dz'  +  E(0).  (5.23) 

4  Jo 

Under  conditions  of  weak  scattering, 

E0(z)**E0  (5.24) 

exponential  attenuation  of  the  sound  in  the  gas, 

5(c)  =  50e-Q*  (5.25) 

and  no  contribution  to  the  scattered  field  for  z  less  than  zero, 

£(0)  =  0  (5.26) 

the  scattered  electric  field  strength  is  given  by 

E(z,  l)  =  0(<  -  1  )50£0e;'(wo+w'-u/)t  f  e-az'J *k!'dz'.  (5.27) 

4  Jo 
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Figure  5— i.  The  k-vector  matching  condition  required  for  collinear  bragg 
scattering. 


If  *  is  detined  to  be  the  following. 

—  as 

= - 

'  An 

and  the  k-vector  matched  case  is  assumed.  A k  =  0.  the  envelope  of  the  scattered 
electric  field  is  given  by 
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(5.2fl) 


The  ratio  of  the  incident  optical  field  to  the  diffracted  field  is  given  by  7.  and 
the  power  diffraction  efficiency  T  is  just  equal  to  the  square  modulus  of  7. 
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r  =  I7H  =  hi2  = 

Lo 


u-  as 


2g\0 


(3.30) 


'Phis  is  the  case  for  A !:  =  0.  If  the  k-vector  matched  case  is  not  assumed,  then 
tlie  integral  in  equation  (5.27)  must  be  evaluated  for  A/.-  non-zero.  If  the  integral  i> 


evaluated  at  a  position  z  far  from  the  interaction  region  (z  »  l/a),  the  upper  limit 
of  the  integral  can  be  replaced  by  infinity  without  effecting  the  result.  The  integral 


becomes 


- 1-—\ 

Jo  \a  —  jAkJ 


(5.31) 


This  is  a  complex  Lorentzian  response.  The  power  diffraction  efficiency  T  for  the 
case  of  finite  mismatch  Ak  is  given  by 


(f  -  1)S0 


(i+m2)’ 


(5,32) 


where  Ak  —  {u  —  ui0)/^'a  and  w o  is  the  acoustic  frequency  for  which  Ak  is  zero. 
Equation  (5.32)  provides  a  relation  between  the  Bragg  scattering  efficiency  T  and  the 
peak  acoustic  condensation  S0  in  the  gas.  It  also  predicts  a  Lorenzian  dependence 
upon  the  frequency  of  the  acoustic  wave  in  the  gas. 

5.3  Heterodyne  Sensitivity 

Since  the  Bragg  scattered  light  is  upshifted  in  frequency,  detection  is  made 
possible  by  optically  mixing  it  in  a  photodetector  with  the  unshifted  light  reflected 
from  the  surface  of  the  solid.  The  optical  mixing  provides  a  very  sensitive  means 
of  measuring  the  amount  of  Bragg  scattered  light.  The  sensitivity  is  large  because 
the  heterodyne  signal  is  proportional  to  the  product  of  the  signal  (shifted)  and 
reference  fields  (unshifted).  When  the  reference  field  is  large  compared  to  the  signal 
field,  the  heterodyne  sensitivity  can  be  many  orders  of  magnitude  better  than  direct 
detection  of  the  scattered  light.  Furthermore,  the  signal  can  be  optically  mixed  to 
a  frequency  where  there  is  essentially  no  interference  or  background.  To  calculate 
the  sensitivity  of  the  probe,  the  noise  contributed  in  the  detection  and  amplification 
process  must  be  taken  into  account.  In  this  section, the  signal  to  noise  of  the  probe 
will  be  calculated  as  a  function  of  the  Bragg  scattering  efficiency. 

The  sensitivity  of  the  heterodyne  detection  will  be  evaluated  under  the 
following  conditions.  The  incident  probe  beam  of  average  power  Pq  is  reflected 


by  the  sound  in  the  gas  with  power  reflection  coefficient  T,  and  reflected  from 
the  sample  surface  with  power  reflection  coefficient  R.  Both  reflected  components 
propagate  collinearly  back  to  the  photodetector  where  they  generate  a  heterodyne 
signal  at  the  frequency  of  the  acoustic  wave  in  the  gas.  The  optical  heterodyne 
power  is  proportional  to  the  product  of  the  fields  from  these  two  reflections.  The 
peak  optical  heterodyne  power  at  the  acoustic  frequency  (tva)  is  given  by 

Phel  =  2s/RPq^/tF0  =  2  \/R\/rPo.  (5.33) 

Conversion  of  the  optical  heterodyne  power  to  an  RF  signal  takes  place  in  the 
photodetector.  The  peak  heterodyne  current  is  given  by 

ffiet  =  QrPheti  (5.3-1) 

where  or  is  the  responsivity  (amps/watt)  of  the  photodetector.  The  average  signal 
power  into  the  load  resistor  is  thus  given  by 

Psignal  =  etRl  =  RV  P0~  Rt .  (5.35) 

There  are  three  sources  which  contribute  noise  to  the  signal  as  it  is  detected 
and  amplified.  They  are  shot  noise,  Johnson  noise,  and  amplifier  noise.  The  total 
noise  power  contributed  by  these  three  sources  is  given  by  [5. 3, 5. 4] 

Pnoise  =  4/\  TDF  +  2earDRP0Ri,  (5.36) 

where  K  is  the  Boltzman’s  constant,  T  is  the  absolute  temperature,  D  is  the 
bandwidth,  F  is  the  noise  figure  of  the  preamplifier,  e  is  the  electron  charge,  ar  is 
the  detector  responsivity  (A/W),  Ri  is  the  load  resistance,  and  Po  is  the  average 
detected  power.  Thus  the  signal-to-noise  ratio  is  given  by 

S_  =  2ar2RrPo2Rl  .. 

N  4/v  TDF  +  2earBRPoRi 

This  equation  establishes  the  relationship  between  the  signal  to  noise  ratio  and  the 
Brag*,  scattering  efficiency  f,  given  the  other  detection  system  parameters. 


One  interesting  point  about  the  Bragg  interaction  must  be  made.  After  the 
incident  probe  beam  is  reflected  from  the  surface  of  the  solid,  it  sees  another 
traveling  acoustic  wave.  The  Bragg  condition  is  met  once  again  for  back  scattering. 
However,  the  frequency  in  this  case  is  down  shifted  because  the  sound  wave  is 
receeding  from  the  solid.  The  down  shifted  light  propagates  back  to  the  solid, 
where  it  is  reflected,  and  then  propagates  away  with  the  upshifted  light.  The  result 
is  a  double  sidebanded  carrier.  The  phase  relation  between  the  upper  and  lower 
side  band  relative  to  the  carrier  is  identical  to  that  found  when  phase  modulation 
takes  place.  The  distinguishing  feature,  however,  is  that  the  low-er  sideband  has 
gone  through  two  reflections  from  the  surface  of  the  solid.  Since  the  reflectivity 
of  the  solid  will  always  be  less  than  one,  the  amplitude  of  the  lower  sideband  will 
always  be  less  than  the  upper  sideband  by  a  factor  of  R,  the  power  reflectivity  of 
the  solid.  The  heterodyne  signal,  with  this  effect  included,  is  now  proportional  to 
the  difference  between  the  sideband  amplitudes  rather  than  the  amplitude  of  the 
single  sideband.  To  include  this  effect  in  the  signal  to  noise  ratio  calculation,  the 
reflectivity  in  the  numerator  of  equation  (5.37)  must  be  replaced  by  the  following 
expression. 

R-R(l-R)2  (5.38) 

One  of  the  attractive  features  of  this  heterodyne  scheme  is  that  an  inter¬ 
ferometer  is  not  required  to  detect  of  the  signal.  This  is  the  result  of  the  fact 
that  the  sidebands  on  the  carrier  do  not  have  equal  amplitudes.  In  contrast,  tech¬ 
niques  which  rely  on  phase  modulation  of  an  optical  carrier,  such  as  displacement 
interferometry,  must  use  a  reference  local  oscillator  which  is  in  quadrature  with 
the  center  carrier  for  demodulation.  These  techniques  generally  require  a  stabilized 
interferometer. 

There  may  be  circumstances  where  the  use  of  an  interferometer  arrangement 
with  the  Bragg  detection  may  be  desirable.  In  cases  where  the  surfaces  are  rough 
or  highly  absorbing,  there  would  be  advantage  to  using  a  reference  which  did  not 
depend  upon  the  surface  reflectivity  of  the  material  being  imaged.  In  those  cases, 


the  probe  beam  could  be  split,  a  portion  sent  into  the  cell  to  be  Bragg  scattered 
by  the  sound  propagating  in  the  gas,  and  the  rest  sent  to  the  reference  arm  of  the 
interferometer.  The  reflection  from  the  sound  propagating  in  the  gas  could  then  be 
mixed  with  the  reference  beam  and  a  signal  obtained  independent  of  the  reflectivity 
of  the  solid.  Under  these  circumstances,  the  signal  to  noise  ratio  would  be  given 
by  equation  (5.37),  without  the  dependence  upon  the  sample  reflectivity.  In  most 
cases,  however,  the  disadvantages  of  using  a  stabilized  interferometer  outweigh  the 
problems  of  having  a  signal  which  is  dependent  upon  the  reflectivity  of  the  sample 
surface. 

5.4  Thermoacoustic  Generation  of  Sound 

When  the  surface  temperature  of  a  solid  varies  harmonically  with  time, 
a  periodic  temperature  distribution  is  generated  in  the  surrounding  gas.  This 
temperature  distribution  couples  to  an  acoustic  wave  through  the  expansion  of  the 
gas.  If  the  amount  of  sound  generated  by  the  expansion  of  the  gas  is  large  compared 
to  the  sound  generated  by  the  thermally  generated  displacement  of  the  solid,  then 
a  measurement  of  the  sound  in  the  gas  is  equivalent  to  a  measure  of  the  ac  surface 
temperature  of  the  solid. 

A  simple  one  dimensional  analysis  can  be  used  to  determine  under  w’hat 
conditions  the  contribution  of  the  gas  will  be  large  compared  to  the  contribution 
of  the  solid.  The  surface  displacement  of  a  solid  when  periodically  heated  at  its 
surface  can  be  shown  to  be  approximately  equal  to  [5.5] 

Ua  =  at0o\tt,  (5.39) 

where  rva  is  the  linear  coefficient  of  expansion  of  the  solid,  0O  is  the  peak  sinusoidal 
temperature  variation  at  the  surface  of  the  solid,  and  X<,  =  \/(‘2k,/u.’)  is  the  thermal 
wavelength  in  the  solid.  This  surface  displacement  U„  causes  an  equal  displacement 
in  the  gas,  since  the  acoustic  displacement  at  the  interface  must  be  continuous. 
In  addition,  the  thermal  expansion  of  the  gas  generates  a  displacement  in  the  gas 


which  is  given  by 


(5.40) 


Ug  —  Otg^O^tgi 


where  the  same  definitions  as  above  apply  for  the  gas.  Since  the  acoustic  power 
density  is  proportional  to  the  square  of  the  displacement,  the  ratio  squared  of  the 
two  displacements  gives  the  relative  amount  of  power  contributed  by  the  expansion 
of  gas  as  compared  with  the  expansion  of  the  solid. 


tftCg 
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(5.41) 


Table  5-1  below  shows  the  relative  power  contribution  of  argon  gas  at  a  pressure  of 
100  atmospheres  as  compared  with  several  solids. 

TABLE  5-1 


Relative  Source  Strength  of  Acoustic  Power  in  the  Gas 


Silicon 

Fused  Quartz 

Gold 

Aluminum 

Exp.  Coeff.  (o„)  2.5  X  10-8 

5.5  X  10"7 

1.3  X  10~5 

2.3  X  10"5 

DilTusivitv  (k)  1.0  X  10-"1 

8.3  X  10“7 

1.3  X  10~* 

9.6  X  10~5 

Power  ratio  {Pg/P«)  4.3  X  10 2 

1.1  X  108 

1.2  X  10l 

5.33 

*References  [5. 6, 5.7] 


As  can  be  seen  from  the  table,  for  many  materials,  the  amount  of  acoustic  power  in 
the  gas  is  essentially  determined  only  by  the  temperature  fluctuation  at  the  surface 
of  the  solid. 

When  the  optical  absorption  depth  Lopl  is  finite,  the  calculation  shown  above 
must  be  modified.  Under  this  condition,  the  displacement  ratio  is  approximated  by 
the  following  expression. 


kjL  _ 

H"  ) 


(5.42) 


In  silicon,  where  the  optical  absorption  depth  is  approximately  4000  A  for  green 
light,  the  power  contributed  by  the  gas  is  nearly  two  orders  of  magnitude  larger 
than  that  contributed  by  the  solid. 
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Since  the  acoustic  power  in  the  gas  is  essentially  determined  by  the  harmonic 
temperature  at  the  surface  of  the  solid,  a  relationship  can  be  established  between 
this  temperature  and  the  acoustic  condensation  in  the  gas.  Because  the  acoustic 
impedance  of  the  solid  is  much  larger  than  that  of  the  gas,  the  thermoacoustic 
generation  of  sound  in  the  gas  can  be  considered  to  be  under  contrained  boundary 
conditions.  Using  the  formalism  from  Chapter  3,  the  acoustic  condensation 
generated  thermoacoustically  in  the  gas  can  be  shown  to  be 


S0  = 


kaCtgOo 


(5.43) 


where  ka  and  kt  are  the  acoustic  and  thermal  k-vectors  in  the  gas,  ag  is  the 
expansion  coefficient  of  the  gas,  and  Oo  is  the  peak  ac  temperature  rise  at  the  surface 
of  the  solid.  This  provides  a  relationship  between  the  magnitude  of  the  sinusoidal 
surface  temperature  of  the  solid  0<>  and  the  magnitude  of  the  acoustic  condensation 
in  the  gas  S0.  This  result  can  now  be  used  to  relate  the  Bragg  scattering  efficiency, 
derived  in  section  5.2.4,  to  the  harmonic  surface  temperature  of  the  solid. 


5.5  Photothermal  Sensitivity 

Using  the  results  from  the  three  previous  sections,  the  photothermal  sensitivity 
of  the  probe  can  be  calculated.  Combining  equation  (5.30)  from  section  5.2  and 
equation  (5.43)  from  the  previous  section,  the  Bragg  scattering  efficiency  T  under 
phase  matched  conditions  can  be  related  to  the  harmonic  surface  temperature  0o. 
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If  this  relation  is  substituted  into  equation  (5.37)  of  section  5.3, 

5  2or2/?rP02/?/ 

N  ~  4KTBF  +  2earBIiP0Ri 


(5.44) 


(5.45) 


an  expression  is  obtained  for  the  signal  to  noise  ratio  as  a  function  of  the  harmonic 
surface  temperature  of  the  solid. 


This  characterizes  the  photothermal  sensitivity  of  the  probe.  Table  5-2  contains  the 
experimental  conditions  under  which  this  sensitivity  will  be  evaluated. 

TABLE  5-2 


Experimental  Conditions  -  Argon  gas  (100  atm.,  300  °K) 


Detector  responsivity 

<*r 

.33  A/W 

Sample  reflectance 

R 

.5 

Incident  Optical  Power 

Po 

5.0  mW 

Load  Resistor 

Pi 

50  n 

Boltzman  constant 

K 

1.38  X  10-23//°K 

Temperature 

T 

300  °K 

Noise  figure 

F 

1.77  (2.5  db) 

Optical  wavelength 

^0 

632  nm 

Dielectric  constant 

c 

1.055 

Acoustic  atten.  const. 

a 

1.7  X  105  /m 

Acoustic  k-vector  (1.03  GHz) 

ka 

1.98  X  107  /m 

Linear  expansion  coefl-. 

<*<> 

1.11  X  lO"3 

Thermal  k-vector  (1.03  GHz) 

kt 

1.2  X  108  /m 

With  these  values,  the  photothermal  sensitivity  is  given  by 


S  f  6.9  X  107  \  2 

N=(— —  p 


(5.17) 


In  a  1  Hz  bandwidth,  the  minimum  detectable  harmonic  temperature  variation  is 
given  by 

9min  =  12  X  lO-^C  (5.18) 


This  sensitivity  compares  very  favorably  with  the  sensitivities  of  other  techniques, 
as  will  be  shown  in  section  5.7.  Furthermore,  this  temperature  variation  need  not 
occur  over  a  region  larger  than  the  optical  probe  spot  size. 


5.8  Acousto-optic  Coupling  Medium  -  Figure  of  Merit 

In  the  previous  section,  a  relation  was  established  between  the  signal  to  noise 
ratio  and  the  ac  temperature  at  the  solid  surface  given  the  optical,  electrical  and 
physical  properties  of  the  detection  system.  This  relation  was  then  used  to  calculate 
the  photothermal  sensitivity,  with  argon  gas  as  the  acousto-optic  coupling  medium. 
Although  argon  gas  is  very  attractive  as  an  acousto-optic  coupler,  many  other  gases 
and  liquids  could  potentially  be  used  to  couple  the  photothermally  generated  sound 
to  the  probe  beam.  In  this  section,  a  figure  of  merit  for  the  coupling  medium  will  be 
derived  from  equation  (5.-16)  of  the  previous  section  and  various  liquids  and  gases 
will  then  be  evaluated  for  use  in  the  photothermal  laser  probe. 

The  sensitivity  of  the  photothermal  probe  depends  upon  the  characteristics  of 
the  acousto-optic  coupling  medium  through  the  following  parameters:  the  dielectric 
constant  e,  acoustic  attenuation  constant  a,  acoustic  velocity  Va,  linear  expansion 
coefficient  ag,  and  thermal  difTusivity  k.  A  figure  of  merit  can  be  defined  using  these 
parameters  in  such  a  way  as  to  make  possible  the  relative  evaluation  of  different 
acousto-optic  coupling  media.  The  signal  to  noise  ratio  calculated  in  the  previous 
section  is  proportional  to  the  properties  of  the  coupling  medium  in  the  following 


S  (f-1  )2a2gK 

N  “ 


(5.18) 


A  figure  of  merit  is  obtained  if  the  right  side  of  this  equation  is  normalized  by  the 
signal  to  noise  ratio  obtained  with  a  reference  medium:  argon  gas  at  a  pressure  of 


100  atmospheres. 


F.M.  =  (8.9  xioV-^f 


(5.19) 


Below  is  found  a  table  of  several  different  gases  and  liquids  with  their  associated 
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figure  of  merit  as  calculated  using  equation  (5.49). 

TABLE  5-3 


Acousto-optic  Coupling  Media 


Media 

Bragg  Freq.  (GHz) 

Figure  of  Merit 

Argon  gas  (100  atm.) 

1.03 

1.0 

Argon  liquid  (87  °K) 

2.7 

360. 

Nitrogen  liquid  (77  °K) 

2.7 

4.7 

Helium  liquid  (4.2  °K) 

0.57 

150 

Water  (300  °K) 

4.9 

0.12 

Air  (300  °K,  1  atm.) 

— -  ■  1  -"= 

1.0 

10~8 

As  can  be  seen,  there  are  reasons  for  using  cryogenic  liquids  as  the  coupling 
media.  This  may  be  especially  interesting  for  evaluating  crystalline  defects  and 
doping  in  silicon  or  other  semiconductors.  At  room  temperature,  the  thermal 
conductivity  of  silicon  is  almost  independent  of  the  defect  concentration  unless  the 
density  is  greater  than  10l8/cni3  (5.8).  The  reason  for  this  is  that  the  dominant 
phonon  scattering  mechanism  at  room  temperature  is  three  phonon  scattering  or 
umklapp  scattering.  Only  at  lower  temperatures  is  the  scattering  dominated  by 
impurity  scattering.  Theory  predicts  that  with  liquid  argon  (87  °Iv)  as  the  coupling 
liquid,  dopant  concentrations  in  silicon  should  be  measurable  down  to  levels  near 
10,8/cm3.  Furthermore,  the  sensitivity  of  the  probe  is  enhanced  by  more  than  two 
orders  of  magnitude  by  the  use  of  liquid  argon  rather  than  the  high  pressure  gas. 
Working  with  liquid  argon  would  also  require  that  a  higher  pump  beam  modulation 
frequency  be  used  as  is  shown  in  Table  5-3. 

5.7  Comparison  of  Photoacoustic  and  Photothermal  Techniques 


It  is  difficult  to  make  a  direct  comparison  of  many  of  the  available  photothermal 
and  photoacoustic  techniques  due  to  the  wide  range  of  applications,  resolution  re¬ 
quirements  and  frequencies  employed.  There  are  several  high  resolution  techniques, 


however,  that  can  be  evaluated  and  compared  in  a  clear  way  when  performing  a 
specific  function.  The  three  techniques  that  will  be  compared  are  photodisplacement 
interferometry,  photoacoustic  detection  with  an  acoustic  lens,  and  photothermal 
Bragg  detection.  The  comparison  can  be  made  by  calculating  for  each  technique 
the  smallest  harmonic  temperature  variation  at  the  surface  of  a  solid  that  can  be 
detected  within  a  given  bandwidth.  The  evaluation  will  be  made  at  a  frequency 
of  1  GHz,  since  the  photothermal  Bragg  detection  is  locked  at  this  frequency  by 
the  Bragg  resonance.  After  calculation  of  the  individual  sensitivities  of  these  three 
techniques,  the  results  will  be  discussed  and  compared. 

To  determine  the  sensitivity  of  photothermal  displacement  interferometry,  a 
periodic  temperature  at  the  surface  of  a  solid  must  be  related  to  the  thermoacous- 
tically  generated  displacement.  Then  using  published  values  of  displacement  sen¬ 
sitivity,  the  sensitivity  to  periodic  temperature  can  be  calculated. 

If  a  solid  is  periodically  heated  at  its  surface,  the  resultant  harmonic  tempera¬ 
ture  varies  over  a  region  of  approximately  one  thermal  wavelength.  The  resultant 
surface  displacement  is  given  approximately  by  [5.5] 

Ut  =  O0a9\tt,  (5.50) 

where  0o  is  the  peak  harmonic  temperature,  a,  is  the  linear  expansion  coefficient, 
and  \lt  is  the  thermal  wavelength  at  the  frequency  of  the  harmonic  temperature. 
This  provides  the  relation  between  harmonic  temperature  and  surface  displace¬ 
ment.  The  minimum  detectable  surface  displacement  with  interferometry  has  been 
reported  as  3  X  10~*  Awith  a  0.1  second  integration  time  [5.9].  The  displacement 
sensitivity  is  therefore  given  by 

Umin  =  — — —  (5.51) 

y/TTz 

The  combination  of  the  two  previous  relations  provides  a  photothermal  sensitivity 
for  displacement  interferometry. 


For  the  case  of  the  acoustic  lens,  the  surface  displacement  must  be  converted 
to  an  acoustic  power  density  in  the  coupling  liquid.  The  assumption  is  made  that 
the  acoustic  loading  of  the  surface  displacement  by  the  coupling  liquid  is  small. 
This  is  a  good  approximation  because  the  impedance  of  most  liquids  is  very  small 
compared  to  the  impedance  of  solids.  Without  loading,  the  power  density  in  the 
liquid  created  by  the  displacement  of  the  solid  is  given  by 

j  =  \ZtJUl  (5.53) 


where  Zi  is  the  acoustic  impedance  of  the  liquid  and  w  is  the  radian  frequency  of 
the  modulation.  Using  equation  (5.50),  this  power  density  can  be  related  to  the 
surface  temperature  of  the  solid. 

2  =  (5-54) 

The  total  acoustic  power  generated  by  this  periodic  temperature  over  a  given  area 
is 

P  =  ^^j(Area).  (5.55) 

The  photothermal  sensitivity  for  the  acoustic  lens  can  now  be  calculated  using 
the  results  of  section  4.3.  If  the  acoustic  detection  efficiency  of  the  lens  is  near  10-2, 
and  the  noise  is  given  by  10-20/£?  Watts-sec,  the  minimum  detectable  acoustic 
power  at  the  focus  is  approximately  given  by 

10~18 

Pmin  =  — 77 — Watts  -  sec.  (5.56) 

ti  z 


The  photothermal  sensitivity  of  the  acoustic  lens  is  therefore  given  by 


@  min  — 


f2 Pmin  ( 

ZiArea  \  wa 
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(5.57) 


The  expressions  for  the  photothermal  sensitivity  of  interferometry  and  acoustic 
lens  techniques  will  now  be  evaluated  using  material  values  for  silicon  and  a  one 
square  micrometer  area.  Under  these  conditions,  the  minimum  detectable  harmonic 
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temperature  for  the  interferometric  technique  is  given  by 


$min  — 


2.2  X  10~2oC 
s/TTz 


(5.58) 


The  photothermal  sensitivity  of  the  acoustic  lens  is  given  by 

4.1  X  10-4oC 


nin  ‘ — 


sfTTz 


(5.59) 


The  sensitivity  of  Bragg  detection  was  calculated  in  section  5.5.  The  result  was 

1.2  X  10-4oC 


@min  — 


v//n 


(5.60) 


Comparing  these  sensitivities,  it  is  clear  that  for  the  case  of  silicon,  the  Bragg 
detection  compares  favorably  with  the  other  techniques.  Because  the  photothermal 
sensitivity  of  the  interferometric  and  acoustic  lens  techniques  depend  upon  the 
expansion  coefficient  and  thermal  wavelength,  the  relative  sensitivities  will  depend 
upon  the  material  being  examined.  However,  two  points  must  be  made.  First,  it  can 
be  shown  that  the  acoustic  lens  technique  will  always  remain  more  sensitive  than  the 
interferometry.  Second,  it  can  be  shown  that  the  photothermal  sensitivity  of  Bragg 
detection  will  always  be  much  greater  than  that  of  displacement  interferometry  if 
the  pressure  of  the  argon  gas  is  greater  than  120  atmospheres. 


5.8  Experimental  Verification  of  the  Theory 


5.8.1  Experimental  Implementation 

A  schematic  diagram  of  the  photothermal  surface  probe  is  shown  in  Figure 
5-5.  A  pump  laser,  modulated  at  1.03  GHz,  is  focused  by  an  optical  microscope 
objective  (9)  inside  the  gas  cell  to  a  diffraction  limited  spot  on  the  sample  surface 
(10).  The  absorbed  optical  power  photothermally  generates  sound  which  propagates 
away  from  the  solid.  The  probe  laser  is  directed  through  a  quarter  wave  plate 
(7)  into  the  cell  by  a  partial  mirror  (3),  and  a  polarizing  beam  splitter  (G),  where 
it  is  focused  to  the  same  location  as  the  pump  beam.  There,  it  is  reflected  at 
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Figure  5-5.  Experimental  Arrangement.  ( 1 ) . ( *2 )  fast  silicon  photodiodes:  (3) 
beamsplitter:  (-1)  spike  interference  filter:  (5)  dicroic  mirror;  (6) 
polarizing  beamsplitter:  (7)  quarter  wave  plate;  (8)  optical  window: 
(0)  focusing  objective;  (10)  sample;  (11)  electrical  feed  through; 
(12)  gas  input;  (13)  focusing  micrometer;  (11)  electromechanical 
scanner. 


the  sample  surface  as  well  as  by  the  sound  wave  propagating  in  the  gas.  These 
two  reflected  components  are  then  recollimated  by  the  optical  objective,  and  pass 
through  the  quarter  wave  plate  a  second  time.  The  two  passes  through  the  quarter 
wave  plate  cause  a  90  degree  rotation  of  the  polarization.  The  reflected  probe 
therefore  propagates  through  the  polarizing  beam  splitter.  It  is  then  reflected  by 
a  dicroic  mirror  (5)  and  finally  passes  through  a  narrow  band  interference  filter 
(4)  to  a  fast  photodiode  (1).  The  interference  filter  isolates  the  photodiode  from 
the  pump  laser  power.  The  cell  is  pressurized  to  50-100  atmospheres  by  simply 
connecting  it  through  a  high  pressure  regulator  to  a  standard  argon  gas  cylinder. 
The  sample  is  positioned  at  the  focus  by  a  micrometer  (13)  from  outside  the  gas 
cell.  An  electro-mechanical  scanner  [5.10]  (14)  inside  the  gas  cell  makes  possible  two 
dimensional  scanning  of  the  sample  with  respect  to  the  pump  and  probe  beams.  A 
reference  signal  at  the  frequency  of  the  axial  mode  spacing  of  the  probe  laser  (wr) 
is  detected  in  a  second  photodetector  (2). 

To  avoid  the  problems  of  RF  interference  that  occur  when  detecting  at  the 
modulation  frequency  of  the  pump  laser  (a/a),  the  two  longitudinal  modes  of  the 
Ile-N'e  probe  laser  are  used  to  optically  mix  the  heterodyne  signal  to  a  frequency 
other  than  that  of  the  acoustic  wave.  This  is  accomplished  by  optically  mixing  the 
Bragg  shifted  signal  from  one  axial  mode  (w0  +  wa)  with  the  unshifted  second  axial 
mode  {u0  +  ojr).  The  result  is  a  heterodyne  signal  out  of  the  optical  detector  at 
the  difference  between  the  acoustic  frequency  and  the  axial  mode  spacing  (wa  —  wr). 
This  signal  is  then  mixed  down  to  an  IF  frequency,  filtered  and  square  law  detected. 
See  Figure  5-6.  For  phase  sensitivity,  the  signal  can  also  be  mixed  down  to  dc  and 
detected  with  a  lock-in  amplifier. 

5.8.2  Frequency  Dependence 

The  frequency  dependence  of  the  Bragg  interaction  was  measured  in  the 
following  way.  A  small  portion  of  the  pump  beam  was  split  off  and  the  modulated 
optical  power  was  monitored  using  a  fast  photodetector  and  an  KF  spectrum 
analyser.  While  changing  the  modulation  frequency  of  the  pump,  the  signal  from 
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igure  5-6.  IIctcTodvno  detection  svstem. 


the  probe  beam  was  measured.  The  probe  signal  was  then  normalized  by  the 
modulated  pump  power  as  measured  with  the  photodiode  and  spectrum  analyser. 
The  experimental  results  are  shown  with  the  theoretically  expected  Lorentzian 
frequency  dependence  in  Figure  5-7. 

5.8.3  Absolute  Sensitivity 

The  absolute  sensitivity  of  the  probe  to  acoustic  power  was  measured  using 
a  1  GHz  acoustic  flat  as  the  source  of  sound  in  the  gas.  The  acoustic  flat  was 
chosen  to  measure  the  sensitivity  for  two  reasons.  First,  the  flat  provides  a  fairly 
uniform  acoustic  power  density  over  a  relatively  large  area  (104//m2).  This  makes 
the  positioning  of  the  probe  beam  a  relatively  simple  task.  The  second  reason  is 
that  it  is  possible  to  calculate  the  acoustic  power  density  in  the  gas  transmitted 
from  the  flat  by  measuring  the  RF  input  power  and  the  insertion  loss  using  a  pulse 
echo  technique.  Measurement  of  the  insertion  loss  of  the  transducer  and  calculation 
of  the  transmission  of  the  sound  from  the  flat  into  the  gas  determines  the  total 
acoustic  power  in  the  gas.  The  acoustic  spot  size  can  be  calculated  from  the  size 
of  the  transducer.  This  information  provides  a  knowledge  of  the  acoustic  power 
density  in  the  gas.  To  avoid  any  miscalculation  due  to  a  resonance  of  the  flat  itself, 
the  insertion  loss  measurements  were  made  using  acoustic  pulses  which  were  shorter 
than  the  round  trip  time  in  the  flat. 

To  locate  the  region  illuminated  by  the  acoustic  power  on  the  front  surface  of 
the  flat,  the  flat  was  illuminated  by  a  focused  modulated  laser.  The  flat  was  scanned 
with  respect  to  the  laser  spot,  and  an  image  was  obtained  of  the  photoacoustic  sound 
detected  by  the  flat  in  receiving  mode.  A  scanned  optical  reflection  image  was  also 
obtained  of  the  same  region.  Comparison  of  the  two  images  identified  the  region 
where  the  acoustic  illumination  was  to  be  found.  The  acoustic  spot  diameter  was 
approximately  100  ;/m.  The  optical  image  was  then  used  to  locate  this  region  for 
positioning  of  the  probe  beam.  Below  are  found  the  experimental  conditions  and 
the  minimum  detectable  acoustic  intensity  /m,„  in  the  gas. 
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A  calculation  of  the  theoretical  sensitivity  to  acoustic  power  in  the  gas  using  the 
experimental  values  above  indicates  that  the  experimental  sensitivity  is  a  factor  of 
3.5  less  than  that  predicted  by  theory.  This  value  is  not  much  larger  than  the 
combined  experimental  uncertainty  in  the  measurement. 

8.8.4  Dependence  on  Pressure 

Using  the  same  experimental  arrangement  described  in  the  previous  section,  the 
pressure  of  the  gas  was  varied  as  the  measurement  was  made.  The  results  showed 
that  the  signal  was  proportional  to  the  pressure  to  a  power  of  2.5  over  a  range  of 
26  to  66  atmospheres. 

S  °c  />“  (5.61) 

The  theory  predicts  that  the  signal  should  vary  as  the  third  power  of  the  pressure  of 
the  gas  [5.11]. 

The  most  plausible  explanation  of  the  discrepancy  between  the  theory  and 
experiment  is  the  following.  Since  the  probe  beam  is  highly  focused  and  the  sound 
delivered  by  the  acoustic  flat  is  essentially  planar,  the  phase  fronts  of  the  two  fields 
are  most  closely  matched  in  the  region  near  the  focus.  Farther  from  the  focus, 
the  spherical  phase  curvature  on  the  optical  beam  makes  efficient  heterodyning 
impossible  due  to  the  mismatch  in  the  phase  curvature  between  the  scattered  and 
reference  beams.  At  low  pressures,  the  acoustic  attenuation  is  high,  so  the  only 
region  which  contributes  to  the  heterodyne  signal  is  near  the  focus.  As  the  pressure 
is  increased,  the  regions  farther  from  focus  begin  to  contribute.  However,  due  to 


the  mismatch  in  the  spacial  phase  fronts,  they  do  not  contribute  as  much  as  would 
be  predicted  by  theory  which  assumes  perfectly  matched  phase  fronts.  This  could 
explain  why  the  measured  pressure  dependence  goes  to  a  power  of  2.5  rather  than 
3.0. 

5.0  Photothermal  Imaging  of  Implanted  Silicon 

To  demonstrate  the  sensitivity  of  the  probe  to  variations  in  the  thermal 
properties  of  materials,  a  silicon  wafer  was  implanted  with  150  IveV  boron  through 
a  resist  pattern  (5.12).  After  implantation,  the  patterned  resist  was  plasma  etched, 
leaving  only  a  thin  oxide.  This  oxide  was  then  removed,  exposing  a  clean  silicon 
surface.  Thus,  the  silicon  was  left  uniformly  implanted,  except  in  the  regions  where 
the  resist  had  isolated  the  silicon  from  the  boron  implant.  In  these  regions,  the 
silicon  was  intrinsic.  The  implant  density  was  calculated  to  be  1018/cm3  at  the 
surface,  increasing  to  1021/cm3  at  a  depth  of  5000  A. 

The  sample  was  imaged  with  the  photothermal  probe,  and  there  was  found  a 
variation  in  the  thermal  response  of  the  implanted  region  (darker)  relative  to  the 
unimplanted  region  (lighter).  The  photothermal  image  is  shown  in  Figure  5-8(a). 
The  spacing  of  the  Il-like  characters  from  center  to  center  is  approximately  15 
microns.  A  3  cc  variation  in  the  optical  reflection  was  measured  from  region  to 
region  due  to  the  effects  of  the  implantion.  This  optical  variation  is  too  small  to 
account  for  the  photothermal  contrast  seen  in  the  image.  However,  to  eliminate  any 
optical  variation,  300  A  of  titanium  were  deposited  on  the  silicon  surface  and  the 
sample  was  imaged  again.  Figure  5-8(b)  again  shows  the  thermal  contrast  due  to 
the  boron  implant.  The  titanium  covered  silicon  region  (left)  actually  generates  a 
larger  photothermal  signal  than  the  silicon  alone.  This  is  largely  due  to  the  greater 
optical  attenuation  constant  of  the  titanium.  The  signal  from  the  uncoated  silicon 
region  (right)  is  not  apparent  in  the  image  because  the  bias  level  and  gain  of  the 
signal  were  adjusted  to  allow  for  maximum  contrast  in  the  titanium/silicon  region. 
The  source  of  contrast  in  the  titanium  covered  silicon  region  is  the  variation  in  the 
thermal  wave  reflectance  at  the  interface  between  the  titanium  and  silicon.  These 
images  clearly  demonstrate  the  sensitivity  of  the  probe  to  changes  in  the  thermal 
properties  of  the  silicon  created  by  the  boron  implantation. 


Figure  5-8. 


5.10  Summary 


The  theory  behind  the  photothermal  laser  probe  has  been  presented  along 
with  an  experimental  verification  of  the  theory.  The  sensitivity  has  been  shown  to 
compare  favorably  with  other  high  resolution  techniques.  A  figure  of  merit  has  been 
defined  and  several  coupling  media  have  been  evaluated.  High  resolution  images 
have  been  obtained  demonstrating  the  resolution  and  sensitivity  of  the  probe  to 
variations  in  the  thermal  properties  of  materials. 

The  features  which  distinguish  this  photothermal  probe  from  other  techniques 
are  the  high  frequency  of  modulation,  the  non-contacting  nature  of  the  evaluation 
and  the  high  sensitivity  to  ac  temperature  variations  over  very  small  areas.  In 
applications  where  the  highest  sensitivity  and  resolution  is  desired,  the  probe  should 
provide  a  unique  tool  for  materials  evaluation. 

The  possibilities  for  application  of  the  probe  to  useful  materials  characteriza¬ 
tion  problems  have  barely  been  investigated.  Future  effort  is  planned  in  this  area. 
There  are  several  materials  and  systems  for  which  the  probe  appears  well  suited. 
Examples  of  such  materials  are  polycrystalline  materials  with  anisotropic  grains, 
very  thin  transparent  films  on  absorbing  substrates,  and  semiconductor  surfaces. 
Another  area  of  future  interest  is  the  use  of  the  cryogenic  liquids  as  coupling  media 
for  the  evaluation  of  semiconductors  at  low  temperatures.  It  is  expected  that  the 
probe  will  provide  new  and  useful  information  in  all  of  these  applications. 
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Appendix 


$WATFIV 

C  THIS  PROGRAM  CALCULATES  THE  AC  AND  DC  TEMPERATURE  RESPONSE 
C  AT  THE  SURFACE  OF  AN  ISOTROPIC  SEMI-INFINITE  ABSORBING  SOLID 
C  WHEN  ILLUMINATED  BY  AN  INTENSITY  MODULATED  GAUSSIAN.  OPTICAL 
C  BEAM. 

C 

C 

C  DIMENSION  VARIABLES  AND  ARRAYS 
C 

C  DEFINITION  OF  VARIABLES 

C 

C  W  =  1/e  HALF  WIDTH  OF  OPTICAL  FIELD 
C  KAPPA  =  DIFFUSIVITY  OF  SOLID 
C  ALPHA  =  POWER  ABSORPTION  CONSTANT  OF  SOLID 
C  F  =  POWER  MODULATION  FREQUENCY 
C  KTKERM  =  THERMAL  K  VECTOR 
C  PNOT  =  TOTAL  POWER  IN  OPTICAL  BEAM 
C  INOT  =  INTENSITY  AT  CENTER  OF  GAUSSIAN  BEAM 
C  PPV  =  ABSORBED  OPTICAL  POWER  PER  UNIT  VOLUME 
C  GRIDAB  =  SIZE  OF  THE  RADIAL  GRID  ELEMENT 
C  GRIDC  =  SIZE  OF  THE  Z  GRID  ELEMENT 
C  VE  =  VOLUME  OF  EACH  ELEMENT  IN  SUMMATION 

C  NAB  =  NO.  OF  THERMAL  WAVELENGTHS  SUMMED  OVER  IN  RADIAL  DIRECTION 
C  NC  =  NUMBER  OF  THERMAL  WAVELENGTHS  SUMMED  OVER  IN  Z. 

C  RES  =  RESOLUTION  FACTOR  OF  INTEGRATION  GRID 
C  TEMP(l)  =  SUM  OF  AC  TEMPERATURE  CONTRIBUTIONS 
C  TEMP  (2)  =  SUM  OF  DC  TEMPERATURE  CONTRIBUTIONS 
C  TEMPM1  =  MAGNITUDE  OF  AC  TEMP  VARIATION 
C  TEMPM2  =  MAGNITUDE  OF  DC  TEMP  VARIATION 
C  THCON  =  THERMAL  CONDUCTIVITY  OF  SOLID 
C  RHO  =  DENSITY  OF  SOLID 
C  CAP  =  HEAT  CAPACITY  OF  SOLID 
C  KAPPA  =  THCON/RHO*C 

C  AABB.CC  ARE  DUMMY  INTEGRATION  VARIABLES 

C  A3. C  ARE  ACTUAL  VALUES  OF  POSITION  CORRESPONDING  TO  AA.BB.CC 
C  ABMAX.CMAX  ARE  UPPER  LIMITS  OF  AABB.CC  AS  DETERMINED  BY 
C  THE  GRID  SIZE  AND  REGION  OVER  WHICH  THE  INTEGRATION  WILL  OCCUR. 

C 

C  DIMENSION  VARIABLES 


COMPLEX* 16  TEMP (2) , KTHERM 

INTEGER  AABB , CC . ABMAX . CMAX , RES 

REAL*8  KAPPA. IN0T.LTHERM.GRIDABC2) ,GRIDC(2) 

C 

C  READ  IN  PARAMETERS 

C 

READ  (5,*)  PNOT.W.F, ALPHA 
READ (5.*)  THCON . RHO , CAP 
READ  (5.*)  NAB, NC. RES 
KAPPA  =  THCON/ (RHO*CAP) 

C 

C  OUTPUT  PARAMETERS 

C 

WRITE (6. 10)  PNOT.W.F 

10  FORMAT (025X, ’BEAM  CHARACTERISTICS ’ . //. 

110X, 'PNOT  =’ .E12.3.2X. 'WATTS’./. 

210X.’W=  ’.E12. 3, 2X. ‘METERS’./. 

310X,’F  =  ’.E12.3.2X, ’HERTZ’./) 

WRITE CS . 12)  KAPPA . ALPHA . THCON 
12  FORMAT C025X. ’MATERIAL  CHARACTERISTICS ’ , //. 

110X. 'KAPPA  =  ‘.2X.E12.3./. 

210X. ’ALPHA  =  ’.2X.E12.3./. 

310X. 'THCON  =  ’ . 2X.E12 . 3. //) 

C 

C  SET  GRIDSIZES  -  DEPENDENT  UPON  THE  RELATIVE  SIZE  OF  THE 
C  THERMAL  WAVELENGTH  WITH  RESPECT  TO  THE  SPOT  SIZE  AND  ATTENUATION 
C  LENGTH 
C 

TEMP (1) =CMPLX (0 . , 0 . ) 

TEMP (2) =CM?LX (0 . , 0 . ) 

DO  1111  M=l,2 
IFCM.EQ.2)  F=l.E-9 
RR=FLOAT (RES) 

C1=DSQRT (KAPPA/ (3. 1416*F))/RR 
C2=W/RR 

IF(C1 .LE.C2)  GO  TO  55 
GRIDAB(l)  =  C2 
GO  TO  56 

55  GRIDAB(l)  =  Cl 

56  C3=l .  /  (ALPHA*P.R) 

IF(C3.GT.C1)  GO  TO  57 
GRIDC(l) =C3 


ooo  onn  ono  nan 


GO  TO  58 

57  GRIDC (1)=C1 

58  GRIDAB (2) =W/RR 

GRIDC (2) =1 . / (ALPHA+RR) 

OUTPUT  GRID  DEFINITIONS 


WRITE  (6,11)  GRIDAB (M) , GRIDC (M) . NAB. NC. RES 
11  FORMAT (025X. ’INTEGRATION  VALUES’,//. 

110X. ’GRIDAB  =’,E12.3,2X, ’METERS*,/. 

210X. ’GRIDC  =’.E12.3,2X, ’METERS’,/. 

310X, 'NAB  =  * ,110./. 

410X, 'NC  =  ’,110,/. 

51 OX. 'RESOLUTION  =  ’.110,/) 

CALCULATE  OTHER  VARIABLES 

IN0T=PN0T*2 . 0/ (W**2*3 . 1416) 

KTHERM=DSQRT (3 . 1418*F/KAPPA) * (1 . . 1 . ) 

LTHERM  =  DSQRT (KAPPA/ (3 . 1416*F) ) 

OUTPUT  OTHER  VARIABLES 


WRITE  <6 . 14)  INOT. LTHERM 
14  F0RJAATC25X. ’CALCULATED  VALUES’,//. 

110X. ’INOT  =’ .2X.E12.3,/, 

210X, 'THERMAL  WAVELENGTH  -  LTHERM  =’ .2X.E12.3,//) 


INTEGRATION 

ABMAX= (NAB*RR) +1 
CMAX=  (NC*RR)  +1 
DO  1010  CC=1.CMAX 
C= (FLOAT (CC) - . 5) *GRIDC (M) 

DO  1000  AABB=1.ABMAX 
AB= (FLOAT (AABB) - . 5) +GRIDAB (M) 

CALL  TRAPS (2. 2. 2) 

PPV=ALPHA*IN0T*EXP(-2*(AB**2)/W+*2)*EXP(-ALPHA*C) 
R=SSRT(AB**2+C**2) 

TEMP (M) =TEMP (M) + (PPV/ (R+THCON) ) *CDEXP (-KTHERM*R) * 
1AB+GRIDAB(M) *GRIDC(M) 

1C00  CONTINUE 


OS 
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1010  CONTINUE 

GUTPUT  RESULTS 

WRITE (6. 1030)  ABMAX.CMAX 
1030  FORMAT (1 OX, ’NUMBER  OF  SUMMATIONS’.//. 

110X, 'ABMAX  =  ’ , 2X . 15 , / , 

310X. ’CMAX  =’.2X.I5,//) 

WRITE (6. 1080)  TEMP(M) 

1080  FORMAT (1 OX. ’TEMPERATURE  REAL  =  ' . 2X.E12 . 5. /. 

125X, ’IMAG  =  ’.2X.E12.5,/) 

1111  CONTINUE 

TEMPM1  =  TEMP(l)  *  DCONJG  (TEMP(l)) 

TEMPM1  =  SQRT(TEMPMl) 

TEMPM2  =  TEMP  (2)  +DCONJG  (TEMP  (2) ) 

TEMPM2  =  SQRTCTEMPM2) 

WRITE (6. 1040)  TEMPM1 

1C40  FORMAT (10X. ’AC  TEMP  MAGNITUDE  FLUCTUATION  =  ’.2X.E12.5./) 
WRITE (6, 1041)  7EMPM2 

1C41  FORMAT (10X,  ’DC  TEMP  MAGNITUDE  FLUCTUATION  =  ’  .2X.E12.5,/) 
RATIO  =  TEMPM2/TEMPM1 
WRITE (6. 1042)  RATIO 

1042  FORMAT (10X, ’RATIO  =  ’ , 2X.E12 . 5, /) 

STOP 

END 

C 

C  INPUT  DATA  FGRMAT 
C 

C  PNOT.W.F. ALPHA 
C  THCON.RHO.CAP 
C  NAB. NC. RES 
$DATA 

l.E-3, . 3E-6, 1 .E9.2.44E6 
170. .2320. .712. 


END 

FILMED 
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